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Abstract 

In this article we review the recently discovered asymptotic integrability in the planar 
Af = 4 SYM theory and discuss its breakdown beyond the asymptotic region due to the 
wrapping interactions. We also discuss novel dynamical tests of the AdS/CFT corre- 
spondence one can perform in the special cases when the wrapping interactions may be 
neglected. 
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Chapter 1 
Introduction 



The AdS/CFT correspondence, proposed by Maldacena pQ in 1997, stating that the string 
theory on the AdS§ x S 5 is dual to TV = 4 supersymmetric Yang-Mills in four dimensions, 
has become one of the prime subjects of interest in gauge and string theory. Since it 
is a strong/weak coupling duality, it offers the unique possibility to investigate a four 
dimensional interacting gauge theory beyond perturbation theory. 

In the last few years, thanks to asymptotic integrability, a great progress has been 
made in confirming the above conjecture in the planar limit. This led, at the same time, 
to precious insights and to a deeper understanding of either of the related theories. Very 
briefly the history of the asymptotic integrability in the planar M = 4 SYM theory can be 
summarized as follows. 

In the groundbreaking paper [2J, the one- loop integrability of the dilatation operator in 
certain subsectors of the gauge side of the AdS/CFT correspondence was discovered. Later 
on [3] the complete one-loop dilatation operator has been found and the corresponding one- 
loop Bethe equations were written down [4]. After many non-trivial steps [5], the form 
of the all-loop asymptotic Bethe equations (ABE) was conjectured [6] up to the so called 
dressing factor [7] , which only contributes starting from the four-loop order. Subsequently, 
relying on the crossing equation proposed in [8] and assuming certain transcendentality 
properties, it was possible to uniquely fix this factor [9]. In this way the asymptotic 
spectrum of the planar M = 4 SYM theory has been completely determined. The asym- 
poticity of these equations means that for a generic operator with L constituent fields the 
corresponding anomalous dimension can be calculated correctly up to the 0(g 2L ) order. 

The appearance of the asymptotic integrability suggests that the spectral problem of 
the AdS/CFT correspondence is intimately related to the dynamics of some lower dimen- 
sional model (integrable one-dimensional spin chain? two-dimensional sigma model?). This 
relation is very surprising, in particular, from the perspective that both the string theory 
on AdSc, x S 5 as well as the Af = 4 gauge theory are higher-dimensional quantum theories! 
It is the purpose of this review to explain how this link can be established for "long" oper- 
ators. For operators of arbitrary length such a link cannot be established before revealing 
"true" nature of the lower dimensional model. 

In this article we will lay all the necessary foundations (chapters 1-3) required to present 
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the set of the asymptotic Bethe equations which, up to the wrapping order, determine the 
anomalous dimension of any single trace operator of the planar Af = 4 supersymmetric 
Yang-Mills theory. In chapter 4 we will discuss analytic properties of operators belonging 
to a particular subsector of the theory, the sl(2) sector. The shortest of these operators de- 
scribe the leading breaking of the Bjorken scaling and their leading anomalous dimension at 
unphysical negative values of the spin is determined by the celebrated BFKL equation. In 
chapter 5 we will discuss this relation and use the results of chapter 4 in order to check the 
veracity of the asymptotic Bethe equations beyond the asymptotic region. The asymptotic 
Bethe equations fail this test and therefore must be modified at and beyond the wrapping 
order. On the other hand, since the wrapping order is controlled by the length L, it is 
possible to study non-petrubative behaviour of anomalous dimensions in the limit L — > oo. 
This question is the subject of chapter 6, where novel dynamical tests of the AdS/CFT 
correspondence are presented that one can perform with help of the asymptotic Bethe 
equations. In chapter 7 we discuss an interesting observation relating the Bethe equations 
of the one of the subsectors of the full theory to the Bethe equations of the Hubbard model. 
This relationship, even though valid only to the first three orders of perturbation theory, 
suggests that a well-defined short range model may be capable to account for the wrapping 
interactions. 

Note added: This work is based on the author's PhD thesis submitted to Humboldt Uni- 
versity, Berlin. 

Note added: After this review has been written, the Y-system for planar AdS/CFT corre- 
spondence [10] and TBA equations [H]-[l3] have been proposed. Both sets of equations 
constitute a significant step forward in non-perturbative formulation of the spectral prob- 
lem of the Af = 4 gauge theory. We will, however, not discuss these recent proposals in 
this article. 
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Chapter 2 



The M = 4 Super Yang-Mills Theory 
and the AdS/CFT Correspondence 

In this chapter we will briefly review the gauge theory side of the AdS/CFT duality and 
subsequently formulate the precise correspondence. 

2.1 Gauge Theory Side 

2.1.1 The M = 4 Super Yang-Mills Theory 

The M = 4 super Yang-Mills theory (SYM) is a maximally supersymmetric gauge theory 
in four dimensions without spin-two fields (gravitons). In what follows we will only review 
the aspects of this theory that are necessary to study the integrable structures appearing 
in the planar limit. For more details and references see |14j . 

The constituents of the M = 4 gauge theory are: six scalar fields <3> m , four fermionic 
fields ^ and the gauge fields A^. Anticipating the transformation rules under gauge 
transformations, we will regard the covariant derivative D^, 

V„ := d„ - igA„ VpW := [X>„, W] = d^W - igA^W + igWA^ (2.1) 

rather then the gauge field as the fundamental field. Here, W stands for a representative 
of the fundamental fields 

W=(ZV,*oa,*£,$,n). (2-2) 

The greek indices belong to the Lorenz algebra so(3, 1) = su(2) x su(2), with 
running from 1 to 4 and a, (3, . . . as well as a, 0, . . . taking values 1, 2. The latin indices 
correspond to the R-symmetry algebra so(6) = su(4). The scalar fields are assumed to 
be in the fundamental representation of so(6) and thus m,n, . . . run from 1 to 6, whereas 
fermions transform in the spinor representation and the corresponding indices a,b, . . . take 
values 1, ... ,4. All fields are assumed to be in the adjoint representation of the gauge 
group U(N), and under the local transformations U(x) G U(N) they transform as 

W) -> U\ W\ {U- l )\ , (A,)) i * U{ (AJ* (U-% - ^g- 1 d,U{ (l/- 1 )) , (2.3) 
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with the repeating indices being summed over. In the following we will suppress the U(N) 
indices. The field strength tensor is defined through 

^> = W" 1 [T>» V v \ = d^Av - d u A^ - ig[A„, A v ] (2.4) 

and transforms canonically under local gauge transformations 

^> >-> UF^U- 1 . (2.5) 

The Lagrangian of the Af = 4 SYM theory is composed of terms of dimension four only 

C YU (x,g) = Tt(^^ u + ^V^ n V^ n -^g 2 [^ m ,^ m ,^ n ] 

In the above formula <r M and a m denote the chiral gamma matrices in four und six dimen- 
sions. 

The equations of motion following from this Lagrangian are invariant under the M = 4 
super Poincare symmetry algebra. The infinitesimal bosonic shifts can be parametrized by 
si*,fi = 0,1,2,3 

S B = s tt P ll , (2.6) 
and under this transformation the set (12.21) transforms as 



5*Vp 


igs J~ U i/j 




5 <3> 






u s * aa 










(2.7) 



Here, we have assumed an adjoint action of the generators on the fields 

5 S W:= [5 S ,W] . (2.8) 
The infinitesimal fermionic translations are parametrized by Grassman variables 

4. = <Ql , (2-9) 
and the corresponding field transformations read 
S W V, = igw a a e a ^% 

M>«a = ligaZai c e aP w^ m ^ n \A<^<s<^ 



n a a/3 W b U H 



8 w n = otolAv^ n . (2.10) 



The transformation laws under the action of Q can be obtained by replacing 



^h^^, w^w (2.11) 

in <^m>. 

Using the explicit form of the parametrizations (12.6j) . (12.91) and with help of the equa- 
tions of motion one can derive from (12.7ft . ( 12. 101) and the respective conjugated transfor- 
mations the following algebraic relations between the generators 

{Q a a ,Q b p} = -2ige af} a%$ m , [P„Q a a } = -ige a/3 a^ , 
{Q« a ,Qp b } = -2ige &$ a2$ m , [P^Q^] = -ige &$ afr% a , (2.12) 

Hereafter, we will use the chiral notation only 

W = (V &p , $ afe , V ab , i>l T aP , . (2.13) 
The relation to the canonical notation is given by 

$ m ~ a b ^ ab . (2.14) 

In a similar way we redefine the translation operator P^ (and subsequently also the con- 
formal boost operator K^) 

P^ofP^. (2.15) 

The advantage of this change of notation is the occurrence of the "fermionic" indices of 
su(2) x su(2) and su(4) only. 



2.1.2 The Conformal and Superconformal Symmtery 

The M = 4 gauge theory posses a further classical symmetry, the conformal symmetry 
|15j . What, however, distinguishes this field theory from other massless field theories is 
the preservation of the conformal symmetry after the quantization. A simple consequence 
of this fact is the vanishing of the beta function 

= <•!=<>. (2.16) 

Thus, the charge does not get renormalized and the momentum-energy tensor remains 
traceless. The only divergent quantities are the wave functions of the fields meaning that 
the scaling dimensions of the operators receive quantum corrections. Indeed, the conformal 
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symmetry severely constraints the form of the two-point correlations functions. While for 
any two local scalar fields O a (x\), Ob(x 2 ) the Poincare symmetry implies that 



(6 A (x 1 )d B (x 2 )) = f AB {x 1 - x 2 ) , (2.17) 

with f(x) being an arbitrary scalar function, the conformal symmetry constraint its form 
to 

(0 A (x 1 )d B (x 2 )) = ° AM mg) ■ (2-18) 
\X\ — x 2 \ 

The quantity A (g) is the scaling dimension of an operator and generically receives quantum 
corrections 

A(g) = A + 1 (g). (2.19) 

Here, Ao denotes the canonical dimension and 7(g) the anomalous part. Standard ar- 
guments of the renormalization theory relate 7(g) of an operator to the corresponding 
wave function renormalization j- log Z. The conformal symmetry also constraints up to a 
constant the form of the three-point function 

(6 A ( Xl ) 6 B ( X2 ) 6 c ( Xi) ) = |xi2|A , 4+AB . Ac| ^f + c ifL, |X3l |A C+ A,-A B • (2-20) 

The correlation functions of four and more fields are not fully determined by the conformal 
symmetry. 

Although the dimension of any operator may be found directly from the correlation 
functions (12.181) and (12.201) it is much more advantageous, as it will become clear later, 
to consider the dilatation operator. The latter is one of the generators of the conformal 
algebra and its eigenvalues are precisely the scaling dimensions 

DO A (x) = A(g)d A (x). (2.21) 

The two spacetime symmetry algebras, the conformal algebra and the supersymmetry 
algebra (together with the so (6) flavor algebra), combine to the super conformal algebra 
psu(2, 2|4). The generators and the structure relations of this algebra are presented and 
discussed in chapter |H It is the unparalleled amount of the symmetries at the quantum 
level that is responsible for many interesting properties of the M = 4 gauge theory. 



2.1.3 The 't Hooft Limit 

In the renowned article of 't Hooft [TB] a novel limit for gauge theories with the U(N) 
gauge groups was proposed, namely iV 1— > 00. In this section we will briefly discuss this 
limit for the Af = 4 SYM theory. 

A common normalization of the actions of the Yang-Mills theories is 

S = -|- f d 4 x£v M (x,l). (2.22! 
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For the purpose of the large N expansion, it is convenient to rescale all fields with g = 

W^gW, (2.23) 

leading to 

S = N j ^£ YM ( x ,g). (2.24) 

All fields in the theory are assumed to be in a hermitian adjoint representation of the gauge 
group U(N) and thus can be represented byNxN hermitian matrices To each upper 
index one assigns an incoming and to each lower index an outgoing arrowo 



j — 



In this notation every propagator is depicted by two parallel lines. For gauge invariant 
operators all indices must be contracted resulting in the "fat" Feynman diagrams (see 
figure |2~7L|) . In this case the contribution of each diagram can be written as [16j 



# N 2 - 29 ' (g 2 )' , (2.25) 

where # is a number, g e denotes the genus of the surface spanned by the diagram and £ 
counts the loop oder. Therefore any physical quantity r\ must admit the following expansion 

OO OO j 

3 = 1 9e=0 

with 7]o being the classical contribution. This applies, in particular, to the scaling dimension 
A(g,N) of the operators. 

An interesting limit, as may be seen from (12.261) . emerges when iV — > oo, g YM — » and 

g 2 = 9 ™J? = const. This limit is called the planar limit, due to the fact that only the 
planar diagrams, that is with g e = (see figure EH]), contribute. 

2.1.4 The Physical Operators 

The basis of local physical operators of the theory is spanned by multiple products of the 
single-trace operators 

6 = Tr(^ 2 . . . Q n ) Tr(fi n+1 . . . Q n+m ) . . . . (2.27) 
Taking into account that the covariant derivative must always act on a field one finds 

Slj G {V k F, V k ^, £> fc $, V k $, V k T) . (2.28) 



1 The reader should recall that the adjoint representation of U(N) can be constructed from the tensor 
product of the fundamental and the antifundamental representation. 
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Figure 2.1: An example of a planar and non-planar Feynman diagram. The grey box 
represents a single-trace operator composed of two elementary fields. 



It can be shown, see e.g. [3], that this set of fields, under the assumption that the Lorenz 
indices may not be contracted, is irreducible with respect to Bianchi identities, equations 
of motion, etc. 

In the planar limit the interactions between operators in different traces in (I2.27P are 
suppressed and it is sufficient to consider the single-trace operators 

Tr^Oa . . . ■ (2.29) 

In this article we will confine ourselves to this case. 



2.2 Formulation of the AdS/CFT Correspondence 

The AdS/CFT correspondence states that the M = 4 SYM theory is dynamically equiva- 
lent to the IIB superstring theory on AdS$ x S 5 if the parameters of the both theories are 
identified as follows 

* = lf' ^7 = 4 ^ (2 ' 30) 
In this formula R denotes the radii of AdS 5 and S 5 and g s stands for the string coupling 
constant. Stated more precisely, to each gauge invariant operator of the M = 4 gauge 
theory there corresponds a dual string state such that all associated physical quantities 
coincide. 

Let | O) be an arbitrary state of the string theory and E(g s , the eigenvalue of this 
state with respect to the string Hamiltonian H defined as an operator conjugated to the 
AdSc, time variable 

H\0) = E(jg.,*)\0). (2.31) 
a' 
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By comparing the string and gauge theory symmetry algebras and their representations, 
one concludes that the string Hamiltonian should correspond to the dilatation operator on 
the gauge theory side. Therefore, there should a exist a state 0(x) such that 



DO(x) = A(g,N)0(x) 



(2.32) 



and 




A(g,N). 



(2.33) 




A(0,oo). 



(2.34) 



According to the above formula scaling dimensions of the planar gauge theory should be 
identified with energies of the free string theory! Unfortunately, the quantization of string 
theory on AdS§ x S 5 product space has not been understood so far and therefore a direct 
verification of (12.331) and (12.341) is not feasible. Despite this fact, it is possible to test this 
equality in certain special cases, some of which we will discuss later. 

The strongest version of the AdS/CFT correspondence states that the duality between 
the both theories holds for arbitrary values of N and g ~ A = (?ym^ i n fl2.3Qj) . In 
particular thus for small gauge groups, as for example SU{2) or SU(3), and for strongly 
coupled string theory. A more modest formulation claims the equivalence in the planar 
limit N —>■ oo, <7ym - ► only. One also cannot exclude that the equivalence of the both 
theories holds only asymptotically so that the O (-^) expansion of the gauge and string 
theory agrees and the discrepancy of the results may be seen only after resummation. 
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Chapter 3 

Foundations of Integrability 



The integrability of a physical system has established itself as an important concept in the 
mechanics and quantum mechanics since it usually leads to deep insights into the dynamics 
of the system. Although it does not necessarily imply solvability, integrability puts severe 
constraints on the physical processes excluding in particular any chaotic behavior. 



3 . 1 Integrability 

In this section we will introduce the concept of classical integrability and discuss its possible 
extensions to the quantum case. 

Let i?(q, p) be a Hamiltonian of a system with N- dimensional phase space. The time 
evolution of the system is then determined through the Hamilton equations 

dH ■ , 

Qj = TT~ J = 1,...,N, 

d Pj 
OH 

= j=l,...,N. (3.1) 

For the Hamilton mechanics the definition of the integrability may be formulated in the 
following way. 

Definition 1. A system is integrable iff there exist N algebraically independent constants 
of motion, that is N functions (Li, . . . , L n ), which satisfy the following two conditions 

i) Vi, j : {Li,Lj} = 0, 

{H,Li} = 0. 

Here, { , } denotes the Poisson bracket. A common way to quantize a system is to 
replace the phase space by a Hilbert space, the variables q, p and the Hamiltonian H by 
the operators 

qi-^q, P^P, H i— > H , (3.2) 
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and the Poisson bracket through the commutator 

{,}-[,]• (3-3) 
Consequently, one could try to formulate the quantum definition of integrability as follows. 

Definition 2. A quantized system is integrable iff there exist N independent operators, 
Li, . . . , Ln, that commute with each other and with the Hamiltonian H . 

Unfortunately, this definition, as we will explain below, is not precise enough; see also 
the discussion in [T7] . 

Let L be one of the constants of motion and let us assume that the spectrum of the 
system is not degenerated. It follows from [H, L] = that H and L have common eigen- 
vectors 

H\* j ) = E j \* j ), L\* j )=7 h \* j ), j = l,...,N, (3.4) 

where N is the dimension of the corresponding Hilbert space. Therefore the operator L 
can be decomposed into the sum of the projectors 

N 

L = J2vjPj, with P, = |^>®(^l- (3.5) 
On the other hand, each projector Pj can be related to the Hamiltonian as follows 

P k = f] E—*!l. (3. 6) 

Putting this formula into the decomposition (13.51) one finds 

N N 

L = Y, H J - 1 F H ■■■, E n)m (3-7) 

j=l 1=1 

Any two commuting operators are thus algebraically dependent and the set of operators 
{H, L\, . . . , L n } can be at most linearly independent. This, however, is not much of a 
constraint on the quantum system since for example {H, H 2 , H+H 3 } satisfy this condition. 
There are also models known for which the existence of the additional commuting operators 
Lj was shown without immediate consequences for the understanding of the physics of the 
system. 

A different possibility to define the quantum integrability offer systems which support 
scattering of particles. These need not to be physical, sometimes physical processes can be 
simply interpreted as scattering of certain particles. In what follows, we will discuss such 
one-dimensional systems. 

Let 1 and 2 be two well-separated particles with asymptotic momenta p\ and p^- The 
"incoming" wave function is then given by 

*to(^i,a?2) ~ expi(piXx +p 2 x 2 ) x x < x 2 • (3.8) 
The both momenta satisfy two conservation conditions 
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1. the conservation of energy, e.g. E = ~ (pf + p 2 ) = const, and 

2. the conservation of momentum P = Pi + P2 = const . 

Therefore, after the scattering process the outgoing momenta need to be a permutation of 
the incoming 

Pi = P2, P2 = Pi ■ ( 3 - 9 ) 
The complete asymptotic wave function can be written as 

^ aS ym P {xi,x 2 ) ~ expz(pixi +P2X2) + S lt2 (pi,P2)expi(p 2 x 1 + pix 2 ) x x < x 2 . (3.10) 

The quantity S(pi,p 2 ) defines the corresponding scattering matrix (S-matrix) and can be 
determined using the Hamiltonian. For theories with the underlying translation invariance, 
the S-matrix depends only on the difference of the both momenta S(pi — p 2 ). In the case 
of three and more particles the conservation conditions 1. and 2. are not sufficient to 
determine the state of the system after scattering. Physically, this corresponds to the 
diffractive three- and many-body processes. The n-particle wave function in the asymptotic 
region, 

xi < x 2 < . . . < x n , (3.11) 
can be represented as a sum of the diffractive and the non-diffractive part 

^a Sym p(Xl,X2, . . .,X n ) ~ 2j eXpz(p CT (i)Xl + • • • +Pa{n)X n ) + ^diffractive (^1 , ^2 , • • • , %n) , 

o-en„ 

(3.12) 

where a denotes an element of the permutation group Il n . The non-diffractive part (the 
Bethe ansatz) describes a sequence of consecutive two-body processes. Since the wave 
function for every two-body process can be represented as in ( 13.101) . the coefficients ^(c) 
and ^(cr') that correspond to two permutations a and a' related by a transposition, 



o = {(7(1), (7(2), . . . , a(i), a(i + 1), . . . , a(n)} (3.13) 

a' = {a(l), (7(2), . . . , a{i + 1), «r(i), . . . , , (3.14) 

must satisfy the following relation 

— — = S^i+itiVW'A^+i)) • (3-15) 



Consequently, the non-diffractive part of the wave function is fully determined by the two- 
body S-matrix. This fact allows for the following definition of the quantum integrability 
(see HZ]). 

Definition 3. A quantum system that supports scattering is integrable iff the scattering of 
the particles is non-diffractive 

Vn : V diffractiveixi, x 2 , . . . , x n ) = . (3.16) 
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In particular, this implies that no particles are being created or annihilated in a scat- 
tering process. 

The vanishing of the diffractive part of the wave function is an indication of existence 
of further conserved quantities (charges). A general three-body process is inelastic. If, 
however, in addition to the conservation laws 1. and 2. there exists a supplementary 
condition, e.g. 

Qz = p\ + pI + ■■■+pI = const, (3.17) 

then the outgoing momenta must be again a permutation of the incoming ones. A generic 
n-body scattering process is non-diffractive if there exist n conserved, algebraically inde- 
pendent and symmetric in all momenta quantities Q r = Q r {pi, ■ ■ ■ ,Pn), r = 1, ...,n. The 
first two of these charges are usually assumed to correspond to the total momentum and the 
total energy of the system. The higher charges Q3, ■ ■ ■ ,Q n are eigenvalues of the generators 
Q3, ■ ■ ■ ,Qn of a hidden symmetry which is not manifest at the level of the Hamiltonian 
or the equations of motion. No general method is known how to construct the generators 
Qr, r > 3 for a generic integrable system. For integrable spin chains, however, which we 
will consider in this article one can accomplish this by constructing the transfer matrix or 
with help of the so called "boost" operator, see [T8] . 

Yet another possibility of defining the quantum integrability provides the Lax pair. Let 
L and A be two N x N matrices, with the elements being functions of the operators of the 
system in question such that 

d ^ = i(AL-LA) (3.18) 

is satisfied. The pair (L, A) is called the Lax pair. From the relation (13 . 18[) it follows that 
the time evolution of L{t) is implemented by a unitary transformation generated by A 

L(t) = U(A,t)L(0)U\A,t). (3.19) 

Upon introducing the determinant, D(X) = det (^L(t) — AlJ, one finds that it is time- 
independent 

D(X) = det (t(t) - A f) = det (lI(A, t) [L(0) - A I] rf(A, t)) 

= det (L(0) - A l) = D(0) . (3.20) 

Calogero has shown [19j that the operator D(X) can be unambiguously defined and that 
the following additional relations are satisfied 

H,D(X)]=0, \d(X),D(X')] = 0. (3.21) 

Moreover, D(X) is according to the definition a polynomial of the iV-th order in A 

N 

D(\) = XJ Qi ■ (3-22) 

j=0 
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Substituting the above formula into (13.211) . one finds the commutation relations 



[#,&■] = fe,4] =0 j,k = l,...,N. (3.23) 

Consequently, the existence of a Lax pair allows to show that the system is integrable in 
the sense of the definition [SJ and it provides a method to construct all the higher charges. 
Unfortunately, no general procedure is known how to determine the corresponding L nd 
A matrices. It is worth mentioning that a similar method was used to show [20J that the 
classical equations of motion of the IIB superstring theory on AdS^ x S 5 are integrable. 
This, in turn, allowed to construct the algebraic curve of the AdS/CFT correspondence 

I2H23I- 

All of the attempts to define the quantum integrability discussed above are not general 
enough to be considered as the final definition. Together, however, they portray plausibly 
what integrability means at the quantum level. Since in this article we will mainly discuss 
spin chains, it is convenient to assume the definition [3] as the criterion of integrability. 



3.2 The S-matrix and the Yang-Baxter Equation 

In this section we will generalize the concept of the S-matrix to the scattering processes 
with different type of particles. We will also formulate consistency conditions that S-matrix 
of an integrable system must satisfy. 

Let us consider an integrable system with n particles that according to their "flavor" 
can be divided into I groups. By flavor we mean a value of some additional charge, e.g. 
spin, electrical charge, etc. Let 1 and 2 be two such particles that carry the flavors fi, f 2 
and the momenta k x , k 2 respectively. The following two scattering channels are possiblqj 

{fi, & K, k' 2 } = {f 2 , A, k 2 , h}, {f{, f>, k[, k' 2 } = {/ ls f 2 , k 2 , h} . (3.24) 

The first scattering channel corresponds to the transition, whereas the second to the re- 
flection of the particles. The definition of the corresponding two-body S-matrix must 
incorporate the additional quantum numbers f x and f 2 

Si2(pi,P2)^S f f l f f *(p uP2 ). (3.25) 

A scattering process of n particles may be decomposed into a sequence of two-body scat- 
terings. Since the latter are non- diffract ive, the momenta of the outgoing particles must be 
a permutation of the incoming ones. Any outgoing configuration of the particles may be, 
however, achieved in many physically distinct ways similarly as the ordered set {2,3,1} 
can be obtained by different sequences of transpositions 

{2, 3, 1} = Z x Z 2 Z x Z 2 {1, 2, 3} = Z 2 Z x {1, 2, 3} = . . . . (3.26) 
1 Here we do not consider flavor-changing processes. 
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Here, Zi denotes a transposition of the two neighboring sites i and i + Using the property 

Zf = I, (3.27) 

the both sequences in (13.261) are equivalent if 

Z\ Z 2 Z\ = Z 2 Z\ Z 2 . (3.28) 

The S-matrix must be consistently defined with respect to the identities (13.271) and (13.281) . 
In the case of (13.271) this implies that 

£ sCl^faPi) S tfM^) = 6*6% , (3.29) 

a.b 

where the sum runs over the two possible intermediate states {a, b} = f 2 } and {a, b} = 
{f2,fi}- Physically, equation f!3.29j) implies that two subsequent scattering processes of 
the same particles are equivalent to no scattering at all, see figure [37T1 The identity (13.281) . 
on the other hand, is reflected in the following cubic relation between S-matrices 

Stf'fapa) Sjfcfaps) S a fl b f2 ( Pl , P2 ) = S f A f *( Pl ,p 2 ) ^(pi.pa) S% b f3 (p 2 ,ps) , (3.30) 

with the indices a, b, c, a, b and c being summed oveJl. The relation (13.301) is the celebrated 
Yang-Baxter equations. It implies that in a non-diffractive scattering process the sequence 
of the two-body scattering processes does not matter, see figure 13.11 

The formulas (13.291) and (I3.30P are sufficient^ to define consistently the scattering of n 
particles. This follows from the fact that the permutation group of n elements, II n , may 
be defined with help of the 2-cycles Zi, i = 1, ...,n — 1 introduced above. Two arbitrary 
Zi and Zj obey 

{Zi Zjf^ = I , (3.31) 

where 

( 1 i = j 

p(i,j)= I 3 \i-j\ = l . (3.32) 
( 2 \i-j\>l 

The first two cases correspond to the identities (13.291) and (13.301) . When \i — j\ > 1, it 
follows from the formula (13.311) that 2-cycles commute, which is trivially satisfied by the 
S-matrices 

•<f. •>,./',.,) S^( Pj ,p j+1 ) = S S l f ^p 3+1 ) . (3.33) 



2 In the presence of fermionic particles both (|3.29p and (|3.30|) must be supplemented with extra minus 
signs. 

3 It should be stressed, however, that these two relations do not imply non-diffractive scattering. In- 
deed, examples are known where the two-body S-matrices satisfy both conditions, but the system is not 
integrable. 
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Pi P2 Pi Pi 

Figure 3.1: A schematic representation of the relations (13.291) and (13.301) . 

3.3 Integrable Spin Chains 

In this section we will apply the previously discussed concepts to the case of spin chains. 
The latter constitute an important subgroup among the integrable models. Moreover, 
the asymptotic integrability in the M = 4 gauge theory discussed in chapter 0] is based 
on identifying the dilation operator with a Hamiltonian of an integrable long-range spin 
chain. 

The Hilbert space of a spin chain of length L is a tensor product of L local quantum 
spaces V 

H = V ®V ® ...®V . (3.34) 

Each single quantum space V is a module of the algebra A @. A spin (particle) is an 
element of this module 

Pa^V } (3.35) 

where A denotes the index with respect to the algebra A. On the Hilbert space 7i one 
defines the Hamiltonian H, 

L-l 

H = ^H i>i+l . (3.36) 

i=i 

The operator Hij+x acts only on the i-th and the (i + l)-th quantum spaces 

... <g> Vj^V (8) . . . . (3.37) 

i i+1 

Additionally, the Hamiltonian must commute with the generators of the symmetry algebra 

[H,A]=0. (3.38) 



In the following we will only consider spin chains for which the the algebra A coincides with the 
symmetry algebra of the Hamiltonian 
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In general one can impose two types of the boundary conditions on the spin chains 

i) open boundary conditions if the last spin L interacts only with preceding (L — l)-th 
spin, and 

ii) closed boundary conditions if the last spin L is assumed to interact with both the 
(L — l)-th spin and the first spin in the chain. 

The second case corresponds to a spin chain on a circle and consequently every local 
quantity 6^ must be related to Qj+L and the domain of summation in (13.36}) is extended 
to include L. In this article we will confine ourselves to periodic boundary conditions for 
which Qj = Qj+L- 

The quantum integrability of a spin chain is understood in the sense of definition [2] or 
El Actually, in this case it is fairly simple to prove the equivalence of the both definitions. 
The operators 

Qr-.H^H r = 3,...,L, (3.39) 

which commute with the Hamiltonian H are non-local and they act on r neighboring sites, 
that is Q r acts on 

. . . ® y ® y ® . . . ® y ® .... (3.40) 

i i+1 i+r—1 

A famous example of an integrable spin chain is the su(2) Heisenberg chain, for which the 
spin takes only two values and the corresponding quantum space is assumed to be 

V = C 2 . (3.41) 

The basis in V is spanned by two states: spin +| and spin — | particles, which we will 
denote by | f ) and | J, ) in what follows. The Hamiltonian can be written in the following 
form 

H 



L 1 

X) z • «Vi - x ) > ( 3 - 42 ) 



4 

i=i 

with the corresponding ground state being composed of L spin-up states 

iJ|t)®...||) = 0. (3.43) 

Excited states are obtained by replacing some +| spins by the — | particles^]. For M 
excitations (magnons) the basis is spanned by 

\xx, . . . , x M ) = I t ) ® • • • ® I T } ®| I ) ® I t > ® • • • ® I t ) ®| I ) ® • • ■ (3-44) 

V v ' N v ' 

X\— 1 X2—XI—1 

1 < x\ < x 2 < ■ ■ ■ < xm < L, xi, . . . , xm = 1, 2, ... . (3.45) 



5 We assume here that the number of the spin-up particles nj is bigger or equal then the number of the 
spin-down particles = L — n^. If > one should take | J. ) (8 . . . | J, ) as the reference vacuum. 



20 



Consequently, each eigenvector \^f) can be decomposed in this basis 

|*)= Yl '${x 1 ,...,Xm)\xx,...,x m ) . (3.46) 

l<xi<. ..<xm<L 

It is here where the integrability plays its profound role and constraints the coefficients 
. . . ,Xm) to be determined purely by the Bethe ansatfi Explicitly, 

M 

W(xi,...,x M )= V ; (^)exp(i^p (T(j) x i ) , (3.47) 
with the scattering of two magnons, in accordance with A3. 15j) . being determined by the 

= S j>j+1 ( Pa(jhP(T(j+1) ) = - 1 + eiPaU)+iPaU+1) _ 2 j PaU) • (3-48) 
It is convenient to introduce the so called rapidities 

u k = Kot^ fc = l,2,...,M, (3.49) 

which allow to re-write the S-matrix elements in the algebraic form 

= - — ; • (3-50) 

The periodic boundary conditions imply that the wave functions must obey 

. . . , x M ) = *(a? 2 , • • • , x M , Xx + L), (3.51) 
leading to the celebrated Bethe equations 

«u + i\ l = fr uj^J^ k = 1 , 2 ,...,M. (3.52) 

In the case of the Heisenberg spin chain, the excitations are all of the same type once the 
reference vacuum has been chosen. For the spin chains with generic underlying symmetry 
algebras, the module V is generally /-dimensional. Let us consider such a generic spin 
chain with K\ particles of type 1, K 2 particles of type 2, etc. so that 

K x >K 2 >...>K h K 1 + K 2 + ... + K l = L. (3.53) 



6 One can also reverse this logic and show that the state (|3.46p together with the ansatz (|3.47|) is under 
the condition (f3~48|) an eigenstate of the Hamiltonian (|3.42|) . Moreover, it can also be proven [24] that all 
2 L states can be parametrized like this. 
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Physical choice corresponds to picking up the most numerous particles (of type 1) as the 
reference vacuu nfl 

|0)i = |1)®...®|1), (3.54) 

N v ' 

L 

and to consider the remaining particles as the excitations on this vacuum. This choice, 
however, breaks the underlying symmetry of the module. The excitations transform under 
the residual symmetry. Under the assumption that the spin chain is integrable an arbitrary 
eigenstate of the Hamiltonian can be represented by a partial state \ A\ . . . A k )\ (with Aj = 
2,3,...,/) and a set of momenta {pi, . . . ,p k } 

\$) = S I \A 1 ...A k ) I , (3.55) 

with 

\A, . . . A k )j = Yl e ipi V^ 2 . . . e^\A u 1,: /,.... i k ) x , (3.56) 

l<il<i2<...<ifc<L 

and 



. . . , A k ;h, . . .,i k ) x = |1) ® . . . ® |1) ®\Ai ) ® |1 ) ® . . . ® |1> ®\A 2 ) ® . . . . (3.57) 

V v ' V v ' 

ii— 1 i'2—ii — l 

The operator S 1 is the many-body scattering matrix of the first level, which due to the 
integrability can be written as a sum over all possible permutations 

S : =^^, (3.58) 
°-£n fe 

with the permutation a corresponding to a given sequence of the two-body processes 

Si = II S h > \ A 'i A ih = E i^fillAiA^ . (3.59) 

It should be stressed that the S-matrix (scattering operator) 5*{ 2 is a (/ — l) 2 x (/ — l) 2 - 
dimensional matrix governing the scattering of two excitations of different flavors. The 
above construction, however, can be repeated after the most numerous excitations (of type 
2) are chosen to be reference vacuum of the spin chain with M = L — K\ sites 

|0) n = |2_^2) r . (3.60) 

M 

The remaining excitations Bi — 3, 4, . . . / on this vacuum scatter with the S-matrix of the 
second level S n , which can be decomposed similarly to (13.581) and (13.591) . Repeating this 
procedure (I — 1) times results in clearing entirely the spin chain from excitations. It should 

7 Hereby we assume that this choice is consistent with respect to the Hamiltonian. 
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be noted that the action of the S-matrix on the first level vacuum is trivial, whereas the 
action on the vacuum of the second level results in an non-trivial scalar phasqj. In the first 
approximation (the first level spin chain) one assumes all of the excitations A4, i = 2, . . . , I 
to scatter with the same phase as the particles of type 2. The error committed in this 
way is partially compensated by considering scattering processes on the vacuum of the 
second level, where the remaining (I — 3) excitations scatter with a phase different from 
the one for the preceding level. The requirement of consistency between the first and 
the second level scattering allows to find explicitly the form of the second level S-matrix. 
Recursive application of this method allows to determine the S-matrices of the k-th level, 
with k = 1, 2, — 1. The periodicity conditions for each of these spin chains lead to 
the nested Bethe equations. This method is known in the literature as the nested Bethe 
ansatz [25] . 

The scattering in generic integrable spin chains with arbitrary values of spin or with 
non-compact symmetry algebras can be uniformly described with the concept of the R- 
matrix. In the following we will only consider systems with the underlying translation 
invariance. Let pa and ps £ V be two particles with the spectral parameter u\ and u 2 
respectively. The R-matrix is defined as an operator that commutes the elements of the 
corresponding modules 

PbM ® PaM = R A§( u i ~ Wa) PaM ® p B {u 2 ) . (3.61) 

Since commuting the modules twice should be equivalent to the action of the identity 
operator, we impose 

REi(u2 - ui) R\%{u x - u 2 ) = 5*5 B D . (3.62) 

Similarly to (I3.30p . the consistency of the scattering processes requires the R-matrix to 
obey the Yang-Baxter equation 

Rj G (u 2 - u 3 ) R%c^i ~ "3) Rab(ui - "2) = Rge(ui - u 2 ) Rihux - m 3 ) R E B F c iu 2 - u 3 ) . 

(3.63) 

Finally, the R-matrices acting on separated modules commute with each other. To a spin 
chain with L lattice sites one can associate the monodromy matrix 

{Te(u)}ii\t::.t = RA d\B L {u) R c ctlt:>) ■ ■ ■ R si» , (3-64) 

which can be intuitively interpreted as a quantity describing an external "ghost" particle 
circulating around the spin chain. The external quantum space V e is also assumed to be 
a module of the symmetry algebra A, though this module must not necessarily coincide 
with the physical one V. The trace of the monodromy matrix with respect to the external 
quantum space 

f(u) = {T e (u)}iilil:i- (3.65) 
8 Here, we exclude the processes 2 + 2 i— > A\ + A2, with A\ 7^ A2 7^ 2. 
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defines the transfer matrix. If the external quantum space coincides with the physical one, 
the transfer matrix becomes a generating matrix of the higher conserved charges 



fW (u) =U expiJ2 u^Qr , (3.66) 

r=2 

where U stands for the translation operator 

U = fV\0) = P li2 P 2 , 3 ...P L -i,L. (3.67) 

Here, P a ^ denotes the permutation operator between the a-th and 6-th vector space. In 
particular the Hamiltonian can be represented as follows 



^^((fWM)- 1 -#«(„)) 



(3.68) 



u=0 



As an example, let us consider the XXX* spin chain with the sl(2) symmetry algebra. The 
modules V can be classified according to the value of the spin s/2. Let |\&) be an eigenvector 
of the Hamiltonian corresponding to M excitations with the rapidities Ui % = 1, . . . , M. 
The state is also an eigenstate of the fundamental transfer matrix (I3.66P 



with 



and 



fV\u)\y) = T( f \u, Ul , u M )\*) , (3.69) 

oo 

T^\u, Ux, . . . , um) — U exp i u r ~ l Q r , (3.70) 



r=2 



M . \s\ M 

Uj + % — 



For the XXX| spin chains T^'(u) is a polynomial of the order L in u. On the other 
hand, the eigenvalue T^>(u, ui, . . . , Um) exhibits poles that depend on the Bethe roots 
Uj,j = 1, . . . , M. Requiring these poles to cancel for arbitrary value of L, one finds the 
Bethe equations 

•LtlFj'L TT *±^H±± k=l,2,...,M. (3.72) 

The case s = 1 corresponds to the Heisenberg spin chain discussed above. 

For the XXX « spin chains it is convenient to introduce beside the fundamental transfer 

matrix also the auxiliary transfer matrix (u) by choosing the external quantum space 
to be V e = C 2 . It was shown in [26] that there exist an operator Q(u) such that 

[Q(u),Q(v)} = [Q(u),f^(v)} = (3.73) 
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and 

{u + i S -) L Q{u + i) + {u- i S -) L Q{u-i) = T (a) (it) Q{u) . (3.74) 

According to the relation (13.731) both T^{v) and Q(u) can be diagonalized simultaneously. 
Let |$) be such eigenstate parametrized by the Bethe roots Ut, % = 1, . . . , M. Equation 
(I3.74p then becomes a functional equation relating the eigenvalue of the auxiliary transfer 
matrix to the eigenvalue of the Baxter operator Q(u), 

{u + i^) L Q{u + i) + {u-i^) L Q{u-i) = T^\u)Q(u). (3.75) 

Here, we have suppressed the explicit dependence of the eigenvalues on {ui}. The equation 
(I3.75P is called in the literature the Baxter equation. Under the following ansatz 

M 

Q(u) = Y[(u - U j) (3.76) 

3=1 

one re-derives the Bethe equations (I3.72p . 
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Chapter 4 



Integrability in the M = 4 Super 
Yang-Mills Theory 

The asymptotic integrability in the planar M = 4 SYM theory allows to compare ob- 
servables on the both side of the correspondence that are not affected by the wrapping 
interactions and therefore paves the way for the novel dynamical tests of the correspon- 
dence. In this chapter we will present the building blocks of the asymptotic integrability 
and discuss the asymptotic spectral equations. 

4.1 The p5u(2,2|4) Super Spin Chain 

4.1.1 Lie Superalgebras and Superconformal Alegbras 

In this section we briefly review essential definitions and results of the Lie superalgebras, 
which we will use in what follows. 

A Lie superalgebra A is an algebra with a Z2 grading equipped with a multiplication 
that obeys the following conditions 

1. skew-symmetry 

Va,beA: (a,b) = -(-l) 9ia)9{b) (b,a), (4.1) 

2. generalized Jacobi identity 

Va,b,ceA: {-iy ia ^ c \a, {b, c)) + (-l)^ b \b, (c, a)) + (-l) 9 ^ b \c, (a, b)) , 

(4.2) 

where g(.) denotes the grade of an element^. The subset 

Aq = {a e A : g(a) = 0} (4.3) 
1 In what follows the multiplication of two fermionic elements will be denoted by { , }. 
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is in view of the conditions (14.11) and (14.21) an ordinary Lie algebra. The fermionc component 
oiA 

A! = {a E A : g(a) = 1} (4.4) 
can be interpreted as a linear representation of Ao since 

(A,A)GA. (4.5) 

Definition 4. A simple Lie superalgebra A for which the representation of its even subal- 
gebra Aq on its odd part Ai is completely reducible is called a classical Lie superalgebra. 

It can be proven that for classical Lie superalgebras the representation of Ao on A\ is 
either 

i) irreducible, or 

ii) is a direct sum of two irreducible representations of Aq. 

This feature allows to define Lie superalgebras of, correspondingly, the first and the second 
kind. 

Definition 5. A classical Lie superalgebra A equipped with a non- degenerate bilinear in- 
variant form is called a basic Lie superalgebra. 

Similarly to the case of the classical semi-simple Lie algebras one can introduce for the 
basic Lie superalgebras the Cartan subalgebra H and the root system {E a ^}, i — 1, . . . , n 
such that 

(Hi, Hj) = 0, 

(Hi, £?±a(j)) = ±Mjj E± a (j) , 

(E a {i), -E-a(j)) = SijHi, (4.6) 

where is the Cartan matrix in the {Hi} basis of the Cartan subalgebra. However, 
in contradistinction to the case of Lie algebras, the choice of the Dynkin diagram is not 
unique. This is due to the fact that there exist a residual freedom in the choice of the odd 
(fermionic) roots. A complete classfication of the fundamental classical Lie superalgebras 
was given in [27] . 

A superconformal algebra is defined to be a Lie superalgebra A for which the even part 
Aq contains the conformal subalgebra so(d, 2) that is spinorially represented on the odd 
component A\. 

4.1.2 The psu(2,2|4) Superconformal Algebra and Its Represen- 
tations 

The su(2,2|4) Lie superalgebra is a superalgebra of the first type, which in the Kac clas- 
sification corresponds to the A(3, 3) Lie superalgebra, see [28J for a pedagogical discussion 
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of this issue. The bosonic (even) component is spaned by su(2, 2) © su(4) © u(l). Taking 
the isomorphism su(2, 2) ~ so(4, 2) into account, one concludes that su(2,2|4) is a super- 
conformal algebra in four dimensions (the spinor representation of the so (4, 2) on the odd 
component of the superalgebra will be explicitly constructed below) . The set of generators 
of su(2, 2, |4) is spaned by 

{Lp, L®, Pa/3, K a p, D; R b ; C | Q a a , Q a a, S", S aa } . (4.7) 

Here, the generators carry "fermionic" indices of su(2) and su(4) respectively and are as- 
sumed to transform canonically. The non- vanishing commutation relations of the su(2, 2|4) 
are given by 

{[D,P a ,];-[D,Kj} = {P a ,;K a ^}, (4.8) 
{ [D, Q a a ] ; [D, S£] ; - [D, Q Aa ] ; - [D, S a ?] } = l -{Q a a ; fi? ; Q &a ;S a ^}, (4.9) 
[K *, P ¥ ] = 8?L a s + 5«L% + 6°sfD , (4.10) 

{QaaiQp} = SIP*?, {^Sg} = 6 a b K?\ 

[S a a ,Pp,\ = S-Q $a , [K°*,Qd = Sfe, (4.11) 



{S:,Q b ,} = 5 b a L a p + 5 a ,R b a + l -5%{D-C) ) (4.12) 
= 5 a b L^-5-R a b + \6^(D + C). (4.13) 

It follows from the above relations that the generator C plays the role of the central 
charge and consequently that su(2,2|4) is reducible. The irreducible part of this algebra, 
psu(2,2|4), can be obtained by considering representations with vanishing central charge. 
It should be stressed, however, that the psu(2,2|4) algebra contrary to the su(2, 2|4) does 
not possess the defining 8x8 matrix representation. 

The above commutation relations are also valid after quantization with all the genera- 
tors except for Up, L°i and R b receiving quantum corrections 

psu(2, 2|4) 3 J : J^J(g). (4.14) 

At weak coupling one expects the following expansion 

oo 

J(g) = J2Jv9 2j - (4-15) 

i=i 

At the one-loop level there is a particularly useful choice of the corresponding Dynkin 
diagram of the su(2, 2|4) superalgebra, see figure I4TT1 The positive and negative roots that 
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correspond to this diagram are given by 

J + e {K a ^S: } S a « } L«(a<[3) } L«(a<P),R a b (a<b)} } (4.16) 

J° e {L%(a = P),ma = p),I%(a = b),D,C}, (4.17) 

J- e {P a ,,Q a a ,Qaa,L$ (a >/3),ZUa>P), Rt(a>b)}. (4.18) 



In the M = 4 SYM theory only non-compact and infinite representations of the super- 
conformal algebra are of physical relevance, which is due to the presence of derivatives of 
fields (each field can be differentiated infinitely many times). Any state in the theory can 
be identified through the following set of Dynkin labels 

{A,Si,s 2 ,qi,p,q 2 ,B,L}, (4.19) 

where A is the eigenvalue of the dilatation operator, the weights [si, s 2 ] classify the spinor 
representations of the Lorentz algebra so(3, 1) = su(2) x su(2) and [q%,p, q 2 ] correspond 
to the flavor algebra so (6). The relation of these labels to the eigenvalues of the Cartan 
algebra elements (14.171) is given by the following formulas 

Sl = Ll-L\, s 2 = Z\-l\, (4.20) 

q 1 = R 2 2 -R{, p = R 3 3 -R 2 2 , q 2 = R\-R\. (4.21) 

The remaining quantities B and L are not related to the weights of the su(2, 2|4) superal- 
gebra and can be seen as eigenvalues of the external automorphisms. The length L counts 
the number of fields in the trace (12.291) and is equal to 1 independently of the value of k 
for each field in (12.281) . The hypercharge B, on the other hand, measures the hyperspin 
of the multiplet (12.281) and is a multiple of |. All physical operators of the theory may 
be classified into highest weight multiplets. Each such multiplet is defined by the highest 
weight state (the primary field) O = \0) 

VJ+ : J + \O} = 0. (4.22) 

Other states of a given multiplet can be obtained by collective action of the lowering 
operators J_ on \0) 

J 1 -J 2 ...J k \0). (4.23) 
An important example of a highest weight multiplet constitute the |-BPS operators 

\Z) L = Tr (Z L ) , Z = $34 • (4.24) 

These highest weight states are additionally annihilated by half of the supersymmetry 
generators and consequently their scaling dimension is protected and does not receive 
quantum corrections, see discussion in section 17.11 

For the purpose of representing fields it is very convenient to use the oscillator realiza- 
tion of the psu(2, 2|4). The oscillators are defined as follows 
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a) bosonic (a a , a£) corresponding to the one copy of the su(2) subalgebra, 

b) bosonic (b a ,bt) corresponding to the second copy of the su(2) subalgebra, 

c) fermionic (c a , c' a ) of the su(4) subalgebra , 

and are assumed to obey the following commutation relations 

[a Q , at] = SI [b A , bp = Sj {c<\ 4} = 61 . (4.25) 
The elements of the psu(2, 2|4) Lie superalgebra for g = can be then represented through 



D = 1 + |at a 7 + fbt b 7 , R a b = cjc a - i^c+c c , 



(4.26) 



C = 1 - -a* a 7 + ^b^b 7 - -c\c c } B = l + ^a+a 7 - -b+b 7 , (4.27) 
2 ^ 2 t 2 2 2 

n« - a tr a O- — h^pt P ■ — atb^ 

Va — d a L 5 V«a — q a ? - r » / 3 — ^a"^) (4 28) 

S« = c ta a , = b d c°, K a ? = a a b^. 

Using formulas (I4.25[) one can easily show that the above realization of the generators obeys 
the commutation relations fl4.8l) - fl4.13p . The set of physical excitations on the vacuum state 
(which is annihilated by all "undaggered" oscillators) is spanned by states for which the 
central charge 

-a+a 7 + -b+b 7 - -« 

2 7 2 7 2 

vanishes. One can easily check that this set coincides with the set of the irreducible fields 



C = 1 - -a+a 7 + -b+b 7 - -ctc c (4.29) 



V h T~ (a f ) fe+2 (b^c^lO), 

~ (e^) k+1 (bt)*^) 1 ^), 
X? fe $ ~ (a^b^c^O), 
pfc^ ~ (a f ) fc (b^ 1 (c^O), 

~ (a f ) fc (b^ 2 (c f ) 4 |0). (4.30) 

Unfortunately, a very disadvantageous feature of this definition is that the vacuum state 
itself does not belong to this class 

C\0) = 1. (4.31) 
One way to overcome this difficulty is to replace the c 3 and c 4 oscillators by d„ , d = 1, 2 

dj = c 4 , d\ = c 3 . (4.32) 
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Figure 4.1: Beauty diagram 



In this notation the Z field is annihilated by all oscillators and can therefore serve as a 
vacuum state. The consequence of such redefinition is the breakdown of the su(4) symmetry 
to su(2) x su(2). On the other hand, the central charge can be rewritten in a more 
transparent form 

C = i(iV b + N A ) -\(N a + N c ) , (4.33) 

where N stands for the counting operators (e.g. N a = a^a 7 ). It follows from (14.331) that 
the excitations on the vacuum state \Z L ) := \Z) ® \Z) can only be created pairwise 

S v ' 

L 

{tf) M {ft) M \Z L ), (4.34) 

with A = {a 1; a 2 , c 1; c 2 } and A = {b 1; b 2 , dj, d 2 }. There are consequently 4 x 4 = 16 
fundamental excitations on each lattice site \Z). Each field in (14.301) is either a fundamental 
excitation or can be represented by a multiple excitation (a composition of the fundamental 
excitations) of a lattice site. 



4.1.3 One-loop Bethe Equations 

It was shown in [4] that the one- loop dilatation operator D 2 (see ( I4.15P ) may be identified 
with the Hamiltonian of an integrable psu(2,2|4) super spin chain. The Bethe Ansatz 
techniques for integrable spin chains with simple Lie algebras as the symmetry algebras 
have been developed in (29] , [30] and subsequently generalized to the case of superalgebras 
[3 lj . This allowed Beisert and Staudacher [4 J to write down the one-loop Bethe equations 
of the planar M = 4 SYM theory. 

As it was explained in section [2 .1.41 in the planar limit it is sufficient to consider only the 
single trace operators. The eigenstates of the one- loop dilatation operator Z) 2 \l/ = 7 2 I '(g)\I/ 
are linear combinations of the basis states (I2.28P 

* = ai Tr(fi x . . . n L ) + a 2 Tr(fi; ...tf L ) + .... (4.35) 

Since the length operator commutes at the one-loop order with the dilatation operator, 
only basis states of equal length need to be taken into account. The complete one- loop 
dilatation operator was found in [3]. Remarkably, it acts only on neighboring fields in each 
trac^l 

L 

D 2 = J2Hi,i+i, (4.36) 
i=i 

2 The explicit form of Hi^ + \ can be found in [3]. 
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which allows to identify single trace operators with states of a closed spin chain 



Tr(fii ...tt L )^\tt 1 ...tt L ):= <g> . . . ® |f2 £ ) . (4.37) 
It follows from the cyclicity of the trace that 

|ftx . . . fi fc fi fc+1 . . . f2 L ) = (_i)(n fc+ i...n i n 1 ...n fe )| 0fc+i ^ _ _ ^ ^ (43g) 

where the prefactor (—l)( n M-i- n £,ni. corresponds to the overall sign of the permutations 
of fermionic fields. The integrability of D 2 was shown in [I] by an explicit construction of 
the psu(2,2|4) R-Matrix. 

Let TC be a Hamiltonian of an integrable super spin chain with sl(K\N) symmetry 
algebra. The diagonalisation of this operator can be performed using the techniques of 
nested Bethe Ansatz, as discussed in section 13.31 For arbitrary values of K and iV and 
arbitrary spin representation it is, however, more convenient to use transfer matrices. The 
transfer matrices for such spin chains were constructed in [29] and [30] . The corresponding 
Bethe equations can be written in the following form [3T] 

u 3 + WkA L rr Uj-ui + {M KjKl 
»J-*VkJ l } l }^ j -u l -iM KjKl ' 

where Mk-Ki denotes the symmetric Cartan matrix and Vk is the spin representation 
vector. There are all together M excitations (with respect to the vacuum state), among 
which Ki of the type 1, K 2 of the type 2, etc. With the help of the Bethe roots Uj one can 
parametrize each conserved quantity: 

1. the momentum (the eigenvalue of the translation operator on the lattice, see (I3.67P ) 

e ,p = e i( P1+ ...+p M ) = JJ u 3 + 2 V «i ; (4 40) 

j=l U 3 ~ 2^ K J 

2. the energy 

M , x 
E = cL ± V ( -. ) and (4.41) 



3. the higher conserved charges 



M 



3=1 



Q r = c r L±Y —■ - — ■ . (4.42) 



u 1 + - 2 vk 3 y- 1 ^-w^y- 1 



The constants c and c r depend on the choice of the Dynkin diagram. Physically they 
correspond to the values of the energy and the conserved higher charges of the chosen 
vacuum state. 
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The psu(2,2|4) superalgebra is a real form of the complex sl(4|4) superalgebra and 
therefore the spectral equations of an integrable spin chain with the psu(2,2|4) as the 
underlying symmetry algebra should also be contained in (14.391) . The highest weight state 
corresponding to the Dynkin diagram presented in the previous section is a ^-BPS field 
and as such is a suitable choice of the vacuum. The set of all excitations on this vacuum 
is a closed sector of the theory. Moreover, every highest weight state can be represented 
in the form (14.341) . which implies that this set coincides with the excitation sector of the 
full planar M = 4 SYM theory. The choice of the vacuum, however, breaks the original 
psu(2,2|4) symmetry. The residual symmetry is 



su(2|2) ®su(2|2) 



(4.43) 



and transforms the sets a 2 , c\, c 2 } and {h\, b 2 , d\, d 2 } in ( 14.341) respectively. The 
central charge H of both superalgebras should be identified with the anomalous dilatation 
operator 5D = D — D in the following way (see [14] ) 



(4.44) 



2 2 



It should be stressed that the overall central charge of the residual symmetry algebra is 
2H. The first level S-matrix can be decomposed into a product of two su(2|2) S-matrices 



'psu(2,2|4)(Pl,P2) = S (p 1 ,p 2 ) ■ ^ S u(2|2)(Pl,P2) ® S SU (2|2) (pi , p 2 ) j , (4.45) 

where So(Pi>P2) is a scalar function. 

The form of the 5 , su (2|2) S-matrices in the asymptotic regime have been derived recently 
in [32] ■ The derivation makes use of the central extension of the su(2|2), where the addi- 
tional central elements are related to some braiding element which modifies the coalgebra 
structure, see (33]. This S-matrix is also invariant under a Yangian [31], see also [35] . 

The one-loop Bethe equations can be found directly form (I4.39p . The Cartan matrix 
corresponding to the Dynkin diagram presented in the previous section reads 



M, 



K h Ki 



( - 2 


+1 














\ 


+1 





-1 

















-1 


+2 


-1 

















-1 


+2 


-1 

















-1 


+2 


-1 

















-1 





+1 


V o 














+1 


-2/ 



(4.46) 



whereas the representation vector takes the following form 

^=(0,0,0,1,0,0,0). 



(4.47) 



33 



Thus, the one-loop spectral equations of the planar Af = 4 SYM are given by 

Ki , K 2 



UAj - 



n 

=1. 
~<i 

n 



■J 



Ki 



Ul,k 
U 2 ,k ~ U X ,j + 



n 



Ui,k - U 2 ,j + 



u l,j + % " Ul,k 



u 2,j 



K 3 

2 TT M 2,A: 



U 2 ,k 


- u u ~ 1 


U3,k 


- U U - | 


U3,k 


- M 2,i + | 


UA,k 


- U 3,j ~ | 


U4,k 


- U 3,j + | 


U 5 ,k 


— M 6,i — 2 



n 

2 i= i 



«3, 



^2,fc - "3,7 + 



U 3 ,k 



UAj 



2=1 

TT M 3,fc ~ M2,j ~ 2 TT 

J= l U 3,fc - U2J + | «3,fc - «3J - 1 ^ U 3 ,fc ~ «4,j + \ 

K 3 



U 3 j + I j-j- u 3:k 



n 

5=1 

^6 
A' 7 



TT M 4,fc — UAj + I T-T U^k — U 5 j — 2 

«5,fc - «5,.j + 1 TT M 5,fc - «4J - | 



n M 6,fc — U 7 j + g TT 
„,._„,._i n 



#5 

n 

2 TT M 6,A; 



M 5,i + « TT M 5,fc 



It 



\ U 5 ,k - U 4 ,j + 2 



3 =1 M 6,fc 
if 7 



"2/| "<<> ~ " : '-J ■ 2 

K 2 



-p-i- U 7}k - U 7 j - i yr U 7 , k - U 6 ,j + 2 



(4.48) 



One notes that the above system of equations is symmetric with respect to the equation for 
the momentum-carrying roots u±. This reflects the decomposition of the S-matrix (14.451) . 

Since the trace is cyclic, the total momentum must be a multiplicity of 2tt and conse- 
quently the Bethe equations must be supplied with the condition 



jp 



K 4 

n 

0=1 



Uaj + 



(4.49) 



The excitation numbers {Ki},i = 1, . . . , 7 appearing in the above equations are uniquely 
determined through the Dynkin labels of a stated 



( Kx \ 

K 3 

K 5 
\K 7 ) 



/|A 
A 
A 
A 
A 
A 
V |A 



\{L-B) 
(L — B) 
\{L-B) 



B) 
B) 
B) 



2 Sl 



p- 

\s 2 



\ 



jQ2 
f?2 



/ 



(4.50) 



3 These relations can be derived from the weights of the corresponding oscillators. 
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Figure 4.2: Beauty diagram with the positive simple roots and the corresponding excitation 
numbers. 

As mentioned before, the anomalous part of the dilatation operator can be identified with 
the Hamiltonian of the psu(2,2|4) super spin chain and thus the energies of the states 
of this spin chain should be proportional to the eigenvalues of the dilatation operator. 
Explicitly, after solving the Bethe equation^ the anomalous dimension can be determined 
from 

Ka , \ 

i= 2 ^^(vh-^) +0{gl) ' (451) 

j=l \ U J ' 2 U 3 2 / 

as follows from (I4.4ip . For the choice of the representation vector made in (I4.47P only 
the excitations (magnons) of type 4 carry the momentum. This allows for the following 
interpretation. A magnon is created by exciting the fourth node of the Dynkin diagram, 
which corresponds to the action of the fourth positive root on the vacuum state Tr (Z L ) , sec 
figure 14.21 Exciting subsequently the neighboring nodes changes the flavor of the magnon, 
but does not create any new "particles" . In a similar manner, one can excite further nodes 
(though always those adjacent to the already excited ones) and as a result the following 
inequalities must be satisfied 

K x < K 2 < K 3 < K 4 > K 5 > K 6 > K 7 . (4.52) 

A detailed study of the psu(2,2|4) representations [36] shows that no adjacent excitation 
numbers, Ki and K i+ i, may be equal. 



4.2 Wrapping Interactions 

The anomalous dimensions at higher orders of perturbation theory can be determined in a 
two different manners. The first method amounts to evaluating the correlation functions 
(I2.18P or (12.201) of the corresponding operators to the desired order. Alternatively, as shown 
in [37] , one can derive the dilatation operator using some further constraints (BM N scaling 
[S5] , integrability, the closure of the symmetry algebra, etc.). For the full planar M = 4 
SYM theory this is a very tedious task already at the two-loop order and this method was 
only applied in some closed subsectors [39]-[ 



4 The roots Uj are in general complex. Since the Bethe equations are invariant under complex conjuga- 
tion, the complex solutions form pairs [uj , u* ) rendering (I4.5ip real. 
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As it was explained in section 12.1. 3[ Feynman diagrams of a U (N) gauge theory consist 
of the "fat" propagators and can be classified according to the genus of the associated sur- 
face. In each genus class one can furthermore distinguish the so called wrapping diagrams. 
In what follows, we will discuss these diagrams in the planar limit. 

Let 0(x) be a single trace operator, which for the sake of simplicity we will assume to 
be built out of two scalar fields 

6 = Tr($i$ 2 ) • (4-53) 
The corresponding two-point function can be determined using the well-know formula 

M ^ lw= aw^ (4 . 54) 

and the Wick theorem. In the above expression O denotes the asymptotic field and C int 
the interaction part of the Lagrangian. In order to define the wrapping interactions it 
is convenient to introduce the so called spectator fields see [32] • The fields ip are 

inserted into the trace 0{x) in the following way 

Tr ($i(x)$ 2 (a:)) >-*• Tr (^{x)$i{x)^{x)§ 2 {x)^{x)} . (4.55) 

Similarly ip are inserted into 0(y). The both auxiliary fields ij) und ip should be contracted 
while calculating (14.541) (it should be noted, however, that ip and ip are free fields) and 
their contraction will be graphically represented by 

tp(x) tp(y) = 

i i 

We exclude contractions for which two lines of the spectator fields cross each other. Also 
two adjacent parallel lines are considered to be equivalent. 

A wrapping diagram is defined to be a Feynman diagram such that for every contraction 
of the spectator fields all the lines of the spectator fields cross the lines of the other fields. 
An example of such diagram is shown in figure 14.31 Because of their topological definition, 
wrapping diagrams can contribute starting from the order 0(g 2L ) only since at lower orders 
at least one contraction line of the spectator fields does not cross any contraction line of the 
other fields. In some special cases, the wrapping interactions may be delayed even beyond 
this order, as for example occurs for operators that can be identified within different closed 
subsectors and with different corresponding lengths. 

As mentioned before, it was possible to determine the dilatation operator to the first 
few orders in some simple closed subsectors [39]-[3l]. Thereby it was assumed that the 
loop order i is smaller than the lenght L of a state on which the dilatation operator acts. 
The reason for this are precisely the wrapping interactions, which are highly non-local and 
cannot be determined through the symmetry algebra. Moreover, the form of the wrapping 
interactions is different for different lengths. 
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Figure 4.3: An example of a planar wrapping diagram that contributes to the two-point 
correlation function. The grey rectangle represents a single trace operator, which is com- 
posed from two fundamental fields. 

4.3 Asymptotic Bethe Ansatz 

At higher-loop orders it is currently unknown with Hamiltonian of which integrable spin 
chain one should identify the dilatation operator. Despite this fact, based on the analysis 
of the closed subsectors [39], [4"0] . a conjecture on the form of the asymptotic all-loop Bethe 
equations for the planar Af = 4 SYM theory was made in [6]. These were subsequently 
proven in [32] under the assumption of factorizability of the S-matrix. 

The Dynkin diagram (called "Beauty") that we have considered so far is, however, 
unsuitable for the purpose of formulating the asymptotic all-loop spectral equations. The 
reason being the fact that the so (6) subsector, which can be consistently obtained by 
truncating all except the three middle nodes of the diagram, is not a closed subsector 
beyond the one-loop level. Let X, y and Z denote the fields of the so(6) subsector. The 
following mixing process 

XyZ* — >WV, (4.56) 

where U, V denote fermionic fields, must be taken into account already at the two-loop 
order (this process also lowers the length of the state, see discussion below). One concludes 
therefore that so(6) Dynkin diagram should not appear as a subdiagram of the psu(2, 2|4) 
at higher-loop order. 

A suitable all-loop choice of the psu(2,2|4) Dynkin diagram is presented in figure |4~41 
Contrary to the Beauty diagram it contains four rather then two simple positive fermionic 
roots. Another pleasant feature of this choice is that the corresponding vacuum state 
remains unchanged 

\Z) L = Tr (Z L ) . (4.57) 

Similarly to the previous case, there are all together sixteen fundamental excitations. We 
present them in the figure I4.4[ which should be understood as follows. The field X corre- 
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sponds to a single excitation of the main (momentum-carrying) node and generally a state 
in the z-th column and j-th row is obtained by exciting once — nodes left to the central 
node and (J — 1) nodes to the right. 

As in the one-loop case, one associates to each single trace operator a state of a closed 
spin chain. The all-loop system exhibits, in contradistinction to the usual spin chains, 
many novel features. First of all, the length of the spin chain must be smaller then the 
loop-order, else the non-local wrapping interactions need to be taken into account and 
their influence on the dynamics of the system is currently unknown. Generically, the £-loop 
Hamiltonian acts simultaneously on £+1 adjacent lattice sites. A feature that rarely occurs 
for integrable systems. In what follows, the regime L > £ will be called the asymptotic 
region and all quantities defined therein will be marked as asymptotic. Probably the most 
striking feature of this asymptotic all-loop spin chain are the fluctuations of length due to 
the flavor mixing processes, e.g. (14.561) . This may seem to invalidate the above definition 
of the asymptotic region, however the supersymmetry delays the wrapping interactions 
for the operators with lower length mixing with the higher length operators. Despite all 
these unusual properties, the asymptotic all-loop spin chain seems to be integrable, though 
this has not been proven rigorously yet. On the other hand, the perturbative asymptotic 
integrability for the first few orders of perturbation theory have been confirmed in some 
subsectors of the full theory, see [43] and [6] . Under the assumption that this holds to all- 
loop order in the asymptotic region, Beisert have derived [32j the corresponding S-matrix 
(I4.45P up to an overall scalar factor. The asymptotic scattering matrix, due to the assumed 
integrability, remains local, in contradistinction to the Hamiltonian, also at higher loop 
orders. Because of the decomposition (I4.45p . it suffices to construct one copy of the su(2|2) 
S-matrix. It turns out, however, that each su(2|2) algebra must be extended with two 
additional local charges in order to overcome the very restrictive particle representations 
of su(2|2). These additional charges vanish on the physical states, which in both cases is 
equivalent to the momentum constraint 1 = e tP = e t<yPl+ '" +PM \ It is interesting to note 
that after introduction of these charges the corresponding S-matrix is uniquely determined 
from the invariance condition 

VJ e su(2|2) k u(l) 2 : [J®I + I®i,^u(2|2)(p 1 ,Pa)] = 0. (4.58) 

As was explained in section I3.3[ the Bethe equations can be obtained from the periodicity 
conditions of the fc-th level spin chain, where in this case k = 1,2,3,4 (the total number 
of excitations is four). 

To formulate these equations it is necessary to introduce, in addition to the rapidity u, 
the deformation variables 



x{u) = l -u |l + yj 1 - j , x ± {u) = x(u ± l -) . (4.59) 
The asymptotic all- loop Bethe equations [6], [32] can then be written as 
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Figure 4.4: All-loop Dynkin diagram of psu(2,2|4) together with the sixteen fundamental 
excitations 
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(4.60) 



The scalar factor a 2 (u,v) appearing on the right hand side of the fourth eqution is the so 
called dressing factor and is closely related to S (pi,p 2 ) in (14.451) 
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x (u)x+(v) 

In [H] it was advocated that its most general form is 

cr 2 (w, v) = exp(2 i 9(u, v)) , 
where the phase 8(u,v) is given through 



(4.6i; 



(4.62) 



9(u k , Uj) = Pr,r+l+2u(g) <lr+l+2v{Uj) - q r {uj) q r +i+2v(u k )) . (4.63) 



r=2 u=0 



In [9], based on the transcendental properties of the scaling function and homogeneity of 
its power expansion in g, the explicit form of the coefficients /3 r>r+ i +2ll (g) was proposed 
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(4.65) 
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C(2/i+l). (4.66) 



It turns out that the dressing factor contributes at weak coupling starting from the four- 
loop order only, see |45j-[47j. The excitation numbers, similarly to the one-loop case, are 
uniquely determined through the labels of the state in question. The explicit relation for 
the all-loop Dynkin diagram reads 
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(4.67) 



The expresion for the anomalous dimension (14.511) is generalized as follows 
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As it was discussed in the preceding section, starting from the order 0(g 2L ) one needs to 
take the wrapping interactions into account and consequently equations (14.601) are gener- 
ically not valid at and beyond this order. In some subsectors, as for example in the case 
of the sl(2) subsector discussed in the next section, the wrapping corrections are delayed 
by supersymmetry and the asymptotic Bethe equations (14.601) remain valid up to the or- 
der O (g 2L+2 ) ■ In this article all quantities that have been calculated with help of the 
asymptotic Bethe ansatz will be marked with the label "ABA" . 
The higher conserved charges are given by 

For symmetric root distribution all odd charges (14.691) vanish. 

Physical solutions of the system (14.601) . that is solutions that correspond to physical 
states of the theory, must have different values of the rapidities at each nesting level 

Vjj^k: u njk ^u nJ n = l,..., 7 (4.70) 

since in the opposite case the corresponding wave function vanishes. 
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Chapter 5 

Analytical Properties of Twist 
Operators 

Twist operators have played a major role in performing tests of the AdS/CFT corre- 
spondence since they are conjectured to be dual to the so called spinning-string on the 
AdSs x S 5 , which in some limit may be treated semiclassically. In this section we will 
review the analytic properties of twist operators. 

5.1 The sl(2) Subsector 

The sl(2) subsector is a closed subsector with operators being composed from scalar fields 
Z and covariant derivatives T> 

Tr (V M Z L ) + ... . (5.1) 

The dots in (15. ip stand for all possible permutations of the derivatives over the Z fields 
with suitable coefficients in front. The length of these operators is equal to the number of 
the Z fields (see subsection 14.1.21) and the covariant derivatives V should be interpreted 
as excitations on the Tr (Z L ) vacuum. The number of excitations is unbounded M = 

0,1,2, and may exceed the length of the operator. This is due to the fact that 

the sl(2) representation is infinite-dimensional. The simplest operators in this sector are 
composed of two scalar fields und arbitrary number of covariant derivatives 

Tr (ZV M Z) +... . (5.2) 

The weights of the s((2) operators can be immediately read off from their field content 
(15. ip . The labels of the primary state must be, however, slightly shifted in order to comply 
with the unitarity [36] 

{A , si, s 2 , q l7 p, q 2 , B, L} = {L + M - 1, M - 1, M - 1, 1, L - 2, 1, 0, L} . (5.3) 

This corresponds to the action of some lowering operators J~ in (I4.18p . A twist of an 
operator is defined as 

T = A - l - (si + s 2 ) . (5.4) 
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According to this definition the twist of these operator is equal L. 

The closure of the st(2) subsector is also reflected in the compact form of the corre- 
sponding Bethe equations. It follows from ( 15.31) together with (14.671) that the excitation 
numbers on the all-loop diagram (figure I4.4D are given by 

{K h K 2 , K 3 , K 4 , K 5 , K e , K 7 } = {0, 0, M - 1, M, M - 1, 0, 0} . (5.5) 

Since the third and the fifth equations can be written as 

P(x 3 , k ) = 0, P(x 5 , k ) = 0, (5.6) 

with 

M M M-l M-l 

p ( x ) = n ( x - x ti) - n ( x - x ij) = c n ^ ~ ^ =c n^ x ^ > ^-^ 

3=1 3=1 3=1 3=1 

one concludes that 

X3,k = X5,k for k = l,...,M-l. (5.8) 
Moreover, it follows from ( 15.71) that 

P( X t,k) _ J~J X 4,k ~ X 4J _ J-J X Xk ~ X 3,A; _ X^ k — £5^ ^ ^ 

P( X i,k) j=l,j^k X ^ k ~ 3=1 3 ' 4 : fc ~ X3 ' k 3=1 a ' 4 . fc ~ a ' 5 ' fc 

This identity allows to reduce the system of three equations to the equation for the main 
roots u^k only 

(+ \ L M — 2 / + — 

On the solutions of this equation one needs to impose the momentum constraint 

M + 

ii : is-") 

fc=i fc 

At the one-loop order equation (15.101) reduces to the Bethe equations of the non-compact 
XXX_i spin chain 

Uk + A (512) 

U k -\) jj^ U k -Uj + l 

The non-compactness is due to the covariant derivatives, which can occur in arbitrary 
number at each lattice site. 
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The Baxter equation, cf. (I3.75p . corresponding to (15.12p takes the following form 

+ Q(« - i) = t(u) Q{u) , (5.13) 

where t[u) = T^ a \u) is the eigenvalue of the auxiliary transfer matrix 

t{u) = 2u L + q L . 2 u L - 2 + ... + q . (5.14) 

The coefficient qL-\ is always zero and g^_ 2 corresponds to the eigenvalue of the Casimir 
opertor 

q L -2 = -Jo(Jo + l)-\L, J = M + h. (5.15) 

The remaining charges q r , r = 0, . . . , L — 3 consitute a complete set of quantum numbers 
that describe uniquely any given state. Equation ( 15.131) . upon fixing the charges, is a dif- 
ference equation of the second order and thus has two algebraically independent solutions. 
One of the solutions is however non-polynomial and we exclude it by writting 

M 

Q(u) = Y[(u - Uj ) . (5.16) 

3=1 

Given a solution, the one-loop anomalous dimension (14.511) may be found from 

7 ABA = 2 9 2 (Q'(~) ~ Q'(~)) + 0(9*) , (5.17) 

as can be easily checked using (15.161) . 

A very interesting feature of twist operators is their scaling behavior for large values of 
the spin. This was first observed in [18] at the one-loop order and was confirmed to hold 
to all orders in (39] and [9]. More precisely, the limit is defined through M — > oo, with 
L growing slower then logarithmically with M. The anomalous dimension to the leading 
order is hence given by 

A-A =i(g) = f(g) logM+... L^oo,M-^oo L<logM. (5.18) 

The universal scaling function f(g) depends only on the coupling constant g and conicides, 
at least up to three- loop order (see |49j), with the L = 2 scaling function. We will discuss 
this limit in detail in chapter [71 



5.2 Twist-Two Operators 

In this section we will discuss a special case of L = 2, in which the anomalous dimension 
may be found in a closed form as a function of M. We start by noting that in this case 
the transfer matrix takes a particularly simple form 

t(u) = 2w 2 -(M 2 + M + i), (5.19) 
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and the difference eqution (I5.13P 

1 

{u + % -f Q{u + i) + (u- % -f Q{u - i) = {2u 2 - (M 2 + M + -)) Q(u) (5.20) 

is solved by the continuos Hahn polynomials 

Q(u) = 3 F 2 {-M, M + 1, i - iu; 1, 1; 1) . (5.21) 

The Bethe roots are simply zeros of this polynomial, see (I5.16p . It is clear that for twist-two 
operators M must take even values since the roots of (I5.2ip do not obey the momentum 
constraint (14.491) for odd M. 

Using the formula (15.171) . one finds 

7 ABA (M) = Sg^M) + O (s 4 ) , (5.22) 

where S\ is the harmonic sum 

M 

( S , 1 (M)=V-. (5.23) 
1 

The formula fl 5 . 2 2 1) determines the anomalous dimension of twist-two operators as a func- 
tion of the spin M; a rare occurrence even at the one-loop level. 

At higher loops it is also possible to reformulate the Bethe equations (15.101) in the form 
of the Baxter equation, see [50J. This allows to expand and subsequently solve the Baxter 
equation order by order in perturbation theory, see [51] . In this article, however, we will 
not present these solutions, but rather concentrate on the corresponding expressions for the 
anomalous dimension. Remarkably, one can guess them in a fairly simple manner using the 
so called transcendentality principle. This method was introduced in [52] and is based on 
the previous observations made in [53] . It assumes that at each order of the perturbation 
theory £ the anomalous dimension is expressed through the generalized hamonic sums of 
the order (2£ — 1), or through the products of zeta functions and harmonic sums for which 
the sum of the arguments of the zeta functions and the orders of the harmonic sums is 
equal to (2£ — 1). The generalized harmonic sums are defined by the following recursive 
procedure [M] 

S . W -£<!gW Sn ^M).±^L Sn „(,•). (5.24) 

The order £ of each sum S ax ,...,a n is given by the sum of the absolute values of its indices 

£ = |di| + . . . \a n \ , (5.25) 

and the order of a product of harmonic sums is equal to the sum of orders of its constituents. 
The canonical basis of the harmonic sums of £-th order is spanned by 



{•San? S a21t a 22 , . . . , S, 



«11 



1 021 1 + I022I = • • • = \an\ + 1 022 1 + • • • + \a u \} , (5.26) 
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where M dependence of the sums is implicit. Each £-ih order product of harmonic sums 
can be decomposed in this basis. 

We will discuss the aforementioned method of determining the higher order corrections 
to the anomalous dimension of twist-two operators by taking the two-loop order as an 
example. In this case, according to the transcendentality principle, the order of the allowed 
harmonic sums is three (the products of zeta functions and harmonic sums do not contribute 
at this order) and thus the basis (I5.26P takes the following form 

f Q O O O O q o q q q 

|D3, D_3, 2 ,1, <J1,2; '-'-2,1; ^1,-2) 2 ,-l, 0_l, 2 , d-2,-1, d-l,-2i 

Si, 1,1, 0-1,1,1, 01,-1,1, Ol,l,-l, 5_i,_i,i, 5"— 1,1, — 1 , Si, -1,-1) 0-1,-1,-1 

} • (5-27) 

It was conjectured in [55] that the index —1 does not appear in the harmonic sums con- 
tributing to the anomalous dimension. Although some Feynman diagrams may lead to 
harmonic sums with index —1, their total contribution should cancel. Moreover, each sum 
for which its first k indices are equal 1 scales in the limit M — ► oo as 

Si,...,i, ak+1 ,..., an (M) ^ #log fc (M) for M » 1 . (5.28) 

One thus concludes that, in view of the scaling properties of the anomalous dimension 
f l5.18p . also 5*11,1 may not contribute to this order. The physical basis of the harmonic 
sums at the two-loop order is consequently a small subset of (15.271) 

{S3, 5-3, 5*2,1, 5i,2, 5_2,i, 5i,_2} • (5.29) 

The two-loop 74 BA anomalous dimension can be expanded in this basis 

7 4 ABA (M) = a 5 3 (M) + c 2 5- 3 (M) + c 3 5 2 ,i(M) + c 4 5 lj2 (M) 

+c 5 5-2,i (M) + c 6 5i,_ 2 (M), (5.30) 

with ci, . . . , Cg being coefficients one still needs to determine. However, it is possible to 
find the two-loop correction to the anomalous dimension for the first few values of M from 
f)5.10p together with (14.681) . This is further simplified due to the fact that the one- loop 
roots are known, see (I5.2ip . Expanding the roots to the two-loop order in (15.101) 

Uk = uf+g 2 u^ + 0{g'), (5.31) 

one finds using (I5.12p a linear system of equations for , which can be easily solved with 
help of any algebra program (e.g. Mathematica). The coefficients Ci, . . . , cq may thus be 
easily determined 

ll^D. = s 3 + 5-3 - 2 fs h2 + 25 2 ,i] - 25i,- 2 . (5.32) 
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It is only a bit more involved to repeat this procedure for the three-loop correction, with 
the result being 



The both expressions f 1 5 . 3 2 j) and (15.331) coincide with the field theory computations, see [55] 
and [53]. The effectiveness of this method can be however appreciated starting from the 
four-loop order, where the usual perturbative calculations become very involved. In [52] 
the four-loop contribution to anomalous dimension was found. We reproduce this result 
in table 15.11 where in the last row the contribution of the dressing phase was printed in 
boldface. In the next chapter we will use this result in order to test the veracity of the 
equations (15.101) since one expects the wrapping interactions to contribute at this order. 

One can also use the above procedure in order to determine the higher conserved 
charges. It turns out that the l-th. loop correction to the 2r-th charge Q2r0 is composed of 
harmonic sums of the (2£ + 2r — 3)-th order. Since harmonic sums of different orders are 
algebraically independent, this confirms the algebraic independence of the higher charges 
and furnishes an indirect proof of the asymptotic integrability. 

5.3 The One-Loop Non-Linear Integral Equation 

In this section we will introduce the concept of the holes which are dual excitations to the 
magnons. With help of the counting function we will rewrite the one-loop Bethe equations 
for the sl(2) subsector in the form of a non-linear integral equation. This new representation 
allows to understand the integrability from the "dual" side and to investigate analytically 
many limiting cases. 

5.3.1 Magnons and Holes 

Using the ansatz (15.161) the left hand side of the Baxter equation (15.131) becomes a poly- 
nomial of the (L + M)-th order and according to the fundamental theorem of algebra it 
must posses (L + M) roots. Putting 




(5.33) 



u = u k 



k 



1,2,. ..,M, 



(5.34) 



All odd charges vanish due to the u 



u symmetry of the root distribution. 
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2.2 

r, 3,l,2,l 



-2 



4 5_ 7 + 6 S 7 + 2 (SL 3l l,3 + 5_ 3 ,2,2 + 5-3,3,1 + 5_ 2 ,4,l) + 3 (S-2,S 

+ 5'-2,3,-2) + 4 (5-2,1,4 — 5-2,-2,-2,1 — 5-2,1,2,-2 — 5-2,2,1,-2 — 5'l,-2,l,3 

— 5l,_2,2,2 — 5*1, -2, 3, l) + 5 (— 5-3,4 + 5-2, -2, -3) + 6 ( — 55,-2 

+ 5^-2,4 — 5_2,-2,l,-2 — 51,-2,-2,-2) + 7 ( — 5_2,-5 + 5_3,_2,-2 

+ 5_2,-3,-2 + 5_2,-2,3) + 8 (5_4,i,2 + 5-4,2,1 — 5_5,_ 2 — 5_4,3 

— 5-2,1,-2,-2 + 5i,_2,l,l,-2) + 9 53,_2,-2 — 10 5i,_2,2,-2 + H 5-3,2,-2 

+ 12 ( — 5_6,1 + 5_2,2,-3 + 5i,4,_2 +54,_2,1 +54,i,_2 —5_3,i,i,_2 —5-2,2,-2,1 

— 5i, 1,2,3 — 5l, 1,3,-2 —5i, 1,3,2 — 5i,2,l,3 — 5i,2,2,-2 — 5l,2,2,2 — 5l, 2,3,1 — Sl,3,l,- 

— 5i,3,i,2 — 5i,3, 2 ,l — 5 2 ,-2,l,2 — 52,-2,2,1 — 52,1,1,3 — 52,1,2,-2 ~ 5 2 , 1,2,2 

Q Q Q Q Q Q Q Q 

~ *->2, 1,3,1 — ^2,2,1,-2 — ^2,2,1,2 _ ^2, 2,2,1 — D 2,3,l,l ~ ^3,1,1,-2 _ ^^,1,1,2 ~ &\ 

— 53,2,l,l) + 13 52,-2,3 — 14 52,-2,1,-2 + 15 (52,3,-2 + 53,2,-2) 
+ 16 (5_4,i,-2 + 5_2,l,-4 —5-2,-2,1,2 — 5-2,-2,2,1 — 5-2,1,-2,2 _ 5_ 2 ,l,l,-3 

— 5i,_3,i,2 — 5i, -3,2,1 — 5i,_ 2 ,-2,2 — 52,-2,-2,1 + 5-2,1,1,-2,1 + "Sl,l,-2,1,- 
+ 5i, i,_2, 1,2 + 5i, i,_2, 2,l) — 17 5_5,2 + 18 ( — 54,-3 — 56,1 + 5i,_3,3) 
+ 20 ( — 5i,_6 — 5i,6 — 54,3 + 5-5,1,1 + 5_4,-2,l + 5_3,_2,2 + 5_2,~ 
+ 5-2,-3,2 + 5i,3,3 + 53,1,3 + 53,3,1 — 5i, 1,-2,3 — 5i,2,-2,-2 — 52,1,-2, 

— 21 £3,4 + 22 (5i,_2,-4 + 52,2,3 + 52,3,2 + 53,-2,2 + 53,2,2) + 23 ( — 5-3,_4 

— 55,2 + 52,-2,-3) + 24 ( — 5_4,_3 + 5i,_4,_2 — 5i,_3,i,_2 — 5i, 1,1,4 — 5i, 1,4,1 

— 5i,3,_2,l — 5i,4,i,i — 53,-2,1,1 — 53,i,_2,l — 54,i,i,i + 5-2,-2,1,1,1 + 5-2,1,-2,1,1 
+ 5^-2,-2,1,1 + 5i, _2, 1,-2,1 + 5i, i,_2, -2,1 + 5i,i,i,_2,-2 + 5i, 1,2, -2,1 + 5i,2, 1,-2,1 
+ 5 2 ,l,l,-2,l) + 25 52,-3,-2 + 26 ( — 52,5 + 5l,4,2 + 52,4,1 + 54,1,2 + 54,2,l) 

+ 28 (5i, 2 ,4 + 5 2 ,1,4 — 5-3,1,-2,1 — 5-2,1,-3,1 — 5i,_2,l,_3) + 30 5_3,i,_3 

+ 32 (5i,5,i + 55,i,i — 5_3,_2,1,1 — 5-2,-3,1,1 — 5i,_3,_2,l — 5i,_2,-3,l 

— 52,2,-2,1 + 5i,2,-2,l,l + 52,1,-2,1,1 _ 5i, 1,1, -2,1, l) + 36 (£1,1,5 + 5l,3,_3 

+ 53,i,_3 — 5i, 1,-3,-2 — 5i, 1,-2,-3 — 5i,i, 2 ,-3 — 5i,2,-2,2 — ^1,2,1,-3 — £2,1,-2,2 

— 52,1,1,-3) +38 5_3,_3,i + 40 ( — 5i, -4,1,1 — 52,-3,1,1 + 5i,i,i,_2,2) 

— 41 £3,-4 + 42 (-£2,-5 + 5i,_ 4 ,2 + 5i,_3,_ 3 ) + 44 (5i,_ 5>1 + 5 2 ,- 3 , 2 + £3,-3,1; 
+ 46 5 2 ,2,-3 + 48 5i,i,_ 3 ,i,i + 60 (5 1 , 1> _ 5 - 5i,i,- 3 , 2 ) + 62 £ 2 ,- 4 ,i + 64 £1,1,1,-3, 

+ 68 (5i,2,-4 + 52,1,-4 — 5i,2,-3,l — 52,l,-3,l) — 72 5i,i,i,_4 — 80 5i, i,_4,i 

— C(3)Si(S 3 -S_ 3 + 2S_ 2 , 1 ). 

ABA 

Table 5.1: The result for see 

256 ' 



-2,-4,1 

-2) 
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the right hand side of (15 . 13[) vanishes and one obtains the Bethe equations (I5.12p . The 
remaining L solutions are roots of the transfer matrix eigenvalue t(u) 



t(u) = 2 \\{u - 4 fe) ) • (5.35) 

k=l 

These additional roots will be called holes in what follows. For u = it^ the Baxter equation 
(I5.13P gives 

W + i\L M (k) _ 

V^J *=l,2,...,i. (5.36) 

The right hand side of this equation contains the product over the usual Bethe roots and 
consequently the hole roots may be found only after the Bethe roots are determined. This 
is related to the fact that the charges qj discussed in the previous section may be expressed 
through the Bethe roots as can be seen from (15.131) together with (15.161) . By comparing 
( I5.14p and (I5.35P one concludes that 

X>h ) = > ( 5 - 37 ) 

3=1 

which corresponds to the momentum constraint (14.491) for the holes. Intuitively, one can 
consider the hole roots as rapidities of the Z fields. We will discuss this point below. 

In what follows we will confine ourselves to the ground states (states with the lowest 
anomalous dimension in the sector) and to even values of the Lorentz spin M, in which 
case both the magnon roots as well as the hole roots are real and symmetrically distributed 
around the origin. We will further assume that L <C M. The charges q k in ( I5.14p depend 
in general on M, however for the assumed values of parameters it was shown in jlS] that 
the second charge q2 is strongly dominating. This implies that two hole roots are much 
bigger then the other 

(!) ( 2 ) fa M , no\ 

In order to derive an integral equation we rewrite the one-loop equations in the logarithmic 
form 

Af-l 
2 

2L arctan (2 Mfc) = 2 7rnfc — 2 arctan (uk — Uj) , (5.39) 

M-l 
J 2 

where we have chosen Arctan branch of the logarithms and grouped the roots as follows 

, 1 3 M-l fKAn . 
u- k = -u k k = — - — . (5.40) 

In [49] it was conjectured that the mode numbers nk of the ground states are given by 

, L-2 .,. , 13 M-l 

n k = k^ — sgn(A;) for k = ±- ± ± — - — . (5.41) 
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The value of the roots grow monotonically with n k . The mode numbers of the hole roots, 
as we will show below using the counting function, can be splitted into two groups. The 
two universal holes have large mode numbers 



Uil L + M-l u o L + M-l 



K 1 = — 2 , = — 2 • (5-42) 

The rest of the hole roots fills the gap around zero opened by the distribution of the magnon 
roots 

nle{-^,...^\. (5.43) 



2 2 

Hence, the following inequlity between magnon and hole roots is satisfied 

\u^\>\u k \>u^ (jVU). (5.44) 

5.3.2 The Counting Function and the Non-Linear Integral Equa- 
tion 

A particularly suitable quantity to describe simultaneously magnons and holes is the scaling 
function, see |56j. In the case of the sl(2) operators it is convenient to adopt the following 
definition (see [57] ) 

M-l 

Z{u) = L(f)(u,-) + (f)(u-u k ,l), where <p(u, f) = i log ( — — — J . (5.45) 

K ~ 2 

The origin of its name is due to the relation 

Z(±oo) = ±tt(L + M) (5.46) 

in conjunction with 

Z( Uj ) = n{2 nj + S-l), J = ±I,...,±^_1, (5.47) 

Z(u^ ] ) = n(2n ( h k) + 5-1), k = l,...,L, (5.48) 

as can be easily confirmed using (15.451) and (I5.12p . Here, we have introduced 

5 = L mod 2 . (5.49) 

Therefore Z(u) is a continuos function which, whenever u is equal to the magnon or the hole 
root, returns the corresponding mode number. With the help of the counting function it is 
straightforward to determine the expression for the mode numbers of the holes. Directly 
from ( 15.411) together with ( 15. 461) - ( T5T481) one confirms the validity of (15.420 and ( 15.431) . 
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The derivation of the non-linear integral equation is based on the following fundamental 
identity (see [56] and references therein) 



k=i j=i 71 

+ r — f( u )\m\og\l + {-l) s e iZ{u+i ^] , (5.50) 

J-oo 7T 

which can be proven by contour integration methods. Applying this identity to the sum 
in (15.451) . one finds after performing Fourier transformation [57] the integral equation for 
the scaling functional 



y< \ ■ r i r(i/2 + » M ) A r 

Z(«) = zLlog— — — — + ^ i \o~ 



r(i/2-i«) ^ r(i(«-«g^ 

+ l im [ a *H.i± log ~ Imlog [1 + (-1)* e ^^°)1 . (5.52) 

a ^°°J- a 71 du T(i(u — v)) 

The non-linear term must be regularised due to the asymptotic behavior of the integrand. 
It should be stressed that the non-linear term may not be integrated by parts, else the 
equation would become linear what contradicts the non-linearity of the Bethe equations. 
The reason for this are discontinuities of the term Imlog [l + (— l) 5 e* z ^ +4 °^] , which need 
to be taken into account while integrating. 

The identity (15.501) allows to express all higher charges (14.691) through the counting 
function Z(u) 

Qp = ~ f ^>)Z(u) - £ + / ^»Imlog [1 + (-1)V^°)] . (5.53) 

J 71 j=1 J 71 

Here, q p (u) denote the corresponding charge densities. The first charge also needs to be 
regularized 

P= lim (- r ^ p '( u )Z(u)-J2p(u^) 

<*-*°° y J_ a ztt ^ 

+ J" ^p'(u) Imlog [l + (-l)V^+ i0 )]J , (5.54) 



2 Duc to superficial divergencies, it is more transparent to apply (|5.50|) to the second derivative of the 
corresponding sum subsequently integrating twice. The integration constants are fixed by antisymmetry 
of Z(u) and 

lim Z'(u) = 0. (5.51) 
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with p(u) = qi(u) being the momentum density of a magnon 

1 it + i/2 

p[u) = - log — . (5.55) 

2 u — i/2 

Using the antisymmetry of Z(u) and p(u) one finds the usual momentum constraint (I4.49[) 

P = 0. (5.56) 
In a similar manner, the one-loop anomalous dimension 72 BA is given by 



7 2 ABA 



4 7e L + 2 ^ {V>(l/2 + iu<i ] ) + ^(1/2 - iu?)} 

+2 £f 4 (*F{i^) toI *[i+<-i>'<''<* M «>] ■ < 5 - 57 > 

where 7e denotes the Euler-Mascheroni constant. 

Equation (I5.52p together with the equation for the holes (15.481) and the mode numbers 
(I5.42p - (l5.43p are completely equivalent to the Bethe equations (I5.12p for the ground states^. 

As it was shown in [57], the last term in (15.521) . which we will subsequently call the 
Imlog-term, at large values of M and for u <C M can be approximated by 

l im F^i* log ^ {U ~ V J> Imlog [1 + (-1) 6 e iZ ^ +i0 A ~ 2 log 2 u . (5.58) 
a^oo J_ a n du Y{i [u — v)) 

Therefore, the equation for the "small" holes (15.481) j = 3, 4, . . . , L takes the following 
explicit form 

The product on the right hand side runs over all hole roots, while the equation itself is valid 
for k = 3, 4, . . . , L. This equation was first derived with the help of the Baxter equation in 
[IB], where also its solutions have been studied in detail. In particular, it was shown that 
the small hole roots (j = 3, 4, . . . L) scale like 

when M > 1. 

All physical quantities can be expressed through hole rapidities and therefore holes can 
be considered as dual excitations of the spin chain. Since the number of holes is equal to the 
number of Z fields, this suggests to identify them with the Z fields. In the dual description 



3 The generalization to all operators in the sector is fairly simple, though we will not consider it in this 
article. 
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the Z fields scatter on the (unphysical) Tr(X> M ) vacuum. Thereby two of the holes move 
very fast and scatter with itself and the remaining slow holes j — 3, 4, . . . , L. It 
is instructive to note that for the twist-two operators there are no slow holes and the 
scattering of the fast holes can be seen as the two-body problem. This clarifies the exact 
solvability in this case. 



5.3.3 The Thermodynamic Limit 

In the case of compact spin chains one usually defines the thermodynamic limit by taking 
the length of the spin chain to be infinite. Non-compact spin chains, in addition, offer 
the possibility to take the number of excitations very large as compared to the length. In 
the context of the AdS/CFT correspondence this is often much more interesting than the 
usual thermodynamic limit. In this section we will assume L — > oo and M — > oo such that 

L<M. (5.61) 

For these values of L and M the mode numbers are explicitly known, see (I5.4ip . 

It was conjectured in [58] that in this limit the magnon roots densly cover the interval 
(— oo, — a) U (a, oo)o Therefore it is appropriate to introduce the density of roots. We 
define it through 

M-l 
1 2 

Po ^ = M ^ Hu-Uj) , (5.62) 

J 2 

from which the normalization condition follows 

dvpo(v) = 1 . (5.63) 



A/-1 
2 



Each sum over the Bethe roots can be expressed through the density 

') 

J 2 

where u(M) denotes the boundary of the magnon distribution 



u(M)\ 

f(u j ) = M[l +1 )dv Po (v)f(v), (5.64) 




lim u(M) = oo . (5.65) 

M— >oo 

Some expressions will need to be regularized with u(M). 

The mode numbers (15.411) can be also obtained from the function^ 

M-l 

M 2 ^ L — 2 
n(u) = — — + 22 0( u ~ u j)^ g — sgn ( M )' n(u k )=n k . (5.66) 

■> 2 



Accordingly, one may show using the properties of orthogonal polynomials that the roots of (|5.21l) 
condense on the whole real axis when M — > oo. 

5 This continuation is, contrary to the counting function, not smooth for finite M. 
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One then notes the following relation between n(u) and the density 



Po(u) = T7-J- ( n ( u ) - ^7r^sgn(w) ) • (5.67) 



M du 

Thus in the limit M -»• oo, cf. CTTI) . 

~Z{u) = 2 7T Po (u) + 2 vr W . (5.68) 

Upon differentiating (I5.45p . one obtains with help of (I5.64p and (I5.68P an integral equation 
for the density 

L — 2 LI ( f~ a ru(M) \ / \ 

2 7rpoH + 2tt— — <5(w) - — — f -2 / + / . m .- ; = 0. 

1 ^ M M U 2 + \ \J- U{ M) J a ) {U ~ Vf + 1 

(5.69) 

This equation must be supplemented with the normalization condition (15.631) . We will 
discuss equation ( 15.691) in detail in section [3 

It follows from (15.601) that in the limit (15.181) three cases should be distinguished. In 
the first one L is kept fixed or goes to infinity slower then logarithmically with respect to 
M; then a = 0. In the second case L = j logM and according to (I5.60p the boundary 
parameter a is constant a = a(j). The last case applies to L diverging faster then loga- 
rithmically, when a is non-constant anymore and grows with M. We will not discuss the 
last case in this review. 



5.4 The Non-Linear Integral Equation at Higher Or- 
ders 

In this section we will generalize the non-linear integral equation to higher orders of the 
perturbation theory. 

Similarly to the one-loop case, the counting function is defined as the logarithm of the 
corresponding Bethe equations (I5.10P 

y, \ ■ t i x ( i / 2 + u ) , • i i + u-u k 
Z(u) = iL\og— — r+l> lo,~ 

nrl i.l ). — ii I ^— ' 



x(i/2—u) i-{u-u k ) 
M 1 j gt M 

~ x(i/2+u)x(i/2-u k ) 



_ 2i Y^ log *(i/2+uMi/2-u fc) +2 J20(u,u k ). (5.70) 



k=l x ^ x(i/2-u)x(i/2+u k ) k=l 



54 



With help of (I5.50P one obtains the desired result 

n ^ , x(i/2 + u) f°° dv ... _ N s 



L poo i 

^2<f>(u-v^\l)- / —(j) , (u-v,l)Im\og[l + (-l) s e iZ ( v+i V] 



dv I . d „ 1 + ^ 



+ / ^ l2^1og- iggy^ -g( tt ,r,) 1 ZQQ 

' ~"~ £c(i/2-u)a; ( i/24-u) 

1 + a! 

+ £ ( 2 .log ^jpZE^ _ M C0) 



DC 



o • " 1 x(i/2+u)x(i/2-v) Q , \ 

2,-log^-j— — -e(«,T,) 

T x(i/2-u)x(i/2+v) 

lmlog[l + (~l) 5 e iZ{v+l0) ] . (5.71) 

Contrary to the previous case (15.521) . linear terms in Z occur on the right hand side of 
this equation. This is due to the fact that the integral kernels are not of a difference form 
anymore and cannot be diagonalized simultaneously with the Fourier transform. 

The mode numbers are not influenced by quantum corrections as long as the coupling 
constant is small and, thus, the relations (15.481) remain valid beyond the one-loop level. 
Equation (I5.7ip together with (I5.48j) is fully equivalent to the asymptotic Bethe equations 
(I5.10P for the ground states. 

The all-loop counting function (15.70!) exhibits similar properties to its one-loop coun- 
terpart. In particular 

^Z(u) = 2 ixp(u) + 2 tt ^^*(«) , (5-72) 

with p(u) being the all-loop density. 

Equation (15.711) can be further simplified by performing a Fourier transformation. 
Defining 

poo 

F(t) = / due- itu F(u), (5.73) 
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one finds after partial diagonalization 

* 2 tt cos (f i^J 2 
w it(e*-l) v y ^ i*(e*-l) e* - 1 w 



+ f£c 0S (^)) 

3=1 7 

- 4/^—/ dt'e->t'K{2gt,2gt')Z(t'), (5.74) 

where £(£) is the Fourier transform of 

£(w) = Imlog [1 + (_i)V z(u+i0) ] . (5.75) 
The kernel K(t,t') is defined through (see [9]) 

t') = + K,(t, + K d (t, t') , (5.76) 

with 

K o{t, t ) = _ , (5.77) 

and 

^^OWMM, (5 78) 

The dressing kernel Kd(t,t') corresponds to the dressing phase (I4.63P and is a convolution 
of the both previous kernels 

K d (t,t') = Sg 2 / dt" •K l (t,2g1?')- p -—k (2gt",t). (5.79) 
Jo e — l 

The asymptotic conserved charges can be expressed through the counting function in a 
similar manner to the one-loop case 



Qp BA = - I ^ A '{v)Z{v) - jr + / ^g p ABA/ (^)Imlog [l + (-l)V^ )] , 

(5.80) 

where g p (w) is the all-loop charge density (14.691) . 

According to the equation (15.741) the function Z(t) exhibits a first-order pole at t — 0. 
The reason for this are the asymptotic properties of the counting function (15.701) . which 
is Fourier transformable only in the principle value sense, similarly to e.g. arctan(u). 
Equation (I5.74p is particularly suitable to investigate different limits in L and M. 
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Chapter 6 



Dynamical Tests of the Asymptotic 
Bethe Ansatz 



In this chapter we will carry out dynamical tests of the asymptotic Bethe equations in the 
sl(2) subsector. For this purpose we will apply the BFKL equation originating from the 
high energy QCD. This equation, which was derived by analyzing leading contributions to 
the hadronic scattering amplitudes, predicts the pole structure of the anomalous dimension 
of twist-two operators. With the help of the four-loop result presented in the previous 
section we will demonstrate that the asymptotic Bethe ansatz is invalid at the wrapping 
order. This supports the hypothesis that if the system remains integrable the wrapping 
interactions must modify the structure of the equations (14.601) . 

6.1 The BFKL Equation and the Double Logarithmic 
Constraints 

An important problem in the theory of hadrons is to determine the behavior of the scat- 
tering amplitudes in the so called Regge limit, when the invariant mass s is much bigger 
then t 



where M denotes the mass scale of hadrons. Regge discovered that in this limit the high- 
energy scattering in quantum mechanics is governed by singularities in angular momentum 
of the partial waves. Later on, this phenomenon was understood to apply to quantum 
field theories as well and in particular to the scattering theory of hadrons. Let A(s,t) be 
a scattering amplitude. The decomposition into partial waves can be performed with help 
of the Mellin transformation, see [59] and [60] for pedagogical reviews, 



s > -t ~ M 2 , 



(6.1) 




(6.2) 



57 



Every pole u of the partial wave A(u,t) contributes to the amplitude A(s,t) the term 



f(t)s"°, (6.3) 

where f(t) is a function of t only. Therefore, knowing all the poles one can reproduce the 
asymptotic behavior of the scattering amplitude. The nearest (dominating) pole is often 
called the Pomeranchuk pole. 

In the limit (16. lft it was shown, see [60J and references therein, that any 2 — > 2 scattering 
in the Regge limit ( 16.11) can be represented in the following form 

A(s,t)=ijr t as ((a log s) m f m , m (t) +a(a log s)^ 1 / m>m _! (*) + ... + a m f mfi (t)\ + O (s°) , 

(6.4) 

where a denotes the coupling constant. The analytic structure in the s plane is related to 
the convergence properties of this series and can be determined only after resummation. 
Generally, for the QCD or the M = 4 SYM, the function f m>n {t) cannot be found explicitly 
and one is forced to make certain approximations. To the leading order one may neglect all 
coefficients except for f m ,m- This is the so called leading logarithmic approximation (LLA). 
In next-to-leading logarithmic approximation (NNLA) one takes additionally f m ,m-i m to 
account, etc. It turns out that in the LLA approximation the hadronic scattering ampli- 
tudes are dominated by interacting gluons propagating in the t channel. These gluons, 
due to their mutual interactions, form collective excitations (see [61]- [64]), the so called 
reggeons. To each reggeon corresponds a sum of infinitely many Feynman diagramms. The 
scattering of two such excitations is equivalent to taking the LLA appoximation and more 
generally one can refine the approximation scheme by taking further reggeons into account. 
It was found in [61]- [64] that the gluon-gluon partial waves obey an integral equation (the 
so called LO BFKL equation), which upon iteration describes the contribution to the 
scattering amplitude coming from a pair of regeeized gluons. Subsequently, the complete 
basis of the homogeneous BFKL equation was deterimined [65J, which allowed to derive a 
functional equation for the eigenvalues. Due to the fact that solely gluons determine the 
leading behavior of the scattering amplitudes, a similar equation may be found in other 
gauge theories. 

Although Af = 4 SYM theory does not have hadrons in its spectrum, it is still possible 
to formulate mathematically the BFKL equation, see [53]. Roughly speaking, the pomeron 
of the Af = 4 gauge theory is described by the non-local gauge-invariant operator 

pomeron = Tr (Z V~ 1+UJ Z ) , (6.5) 

where the parameter uj is assumed to be small. It can be shown [55J that the anomalous 
dimension of this operator (for u — > 0) in the LLA approximation can be found from the 
LO BFKL equation 

*(~Vi +*( 1 + Vi -2tt(l) , (6.6) 



4g 2 \ 2/ V 2 
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where ty(x) is the psi function. Due to the resummation of infinitely many Feynman 
diagramms, the BFKL equation determines non-perturbatively the leading singularities. 
The perturbative solution can be obtained by expanding in 7 

Z^i-gGf (6.7) 

and subsequently substituting the perturbative expansion of the anomalous dimension 

l{g) = 9 2 l2 + # 4 74 + # 6 7e + • • • • (6.8) 
To the first few orders one finds the following perturbative expansion around to = 

It should be noted that the LO BFKL equation predicts the order and the residue of the 
leading pole at each order of the perturbation theory. In order to find the sub-leading poles 
one would need to go beyond the leading approximation scheme. 

The quantum numbers of (16.51) can be formally obtained through analytic continuation 
in M of the corresponding quantum numbers of twist-two operators 

{2,3,...} 3 M M = -1 +u. (6.10) 

In this sense one can also interpret the state (16.51) . This is supported by the well-defined 
analytical properties of the anomalous dimension and higher conserved charges of twist-two 
operators for which closed formulas in M are known. Harmonic sums occurring in these 
expressions, as shown in [66], may be continued analytically to the whole complex plain. 
For example, the one-loop contribution (I5.22j) can be expressed through the psi function 

72 = 8 5!(M) =8(*(M + 1)-*(1)) , (6.11) 

which is defined for MgC. It was shown in [HBJ that all nested harmonic sums (I5.24p can 
be consistently continued to complex values of M by means of psi functions and higher 
transcendental functions. After analytic continuation every harmonic sum of the order I 
exhibits poles at negative integer values of M and the order of the highest pole is always 
smaller or equal £. In particular, the anomalous dimension of twist-two operators is singular 
for M = — 1. This allows to compare the order and the residue of the highest pole with 
the prediction coming from the BFKL equation. We will use this possibility to verify the 
veracity of the four- loop result derived in section 15721 from the asymptotic Bethe equations 

Apart from the BFKL equation, predicting the structure of the poles at M — — 1, there 
also exist constraints following from the double logarithmic behavior of the scattering 
amplitudes, which enable to predict the leading singularities at negative even vaues of M. 
Scattering amplitudes in this limit were studied for QED and QCD in [67], [68] and [69] . 
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According to the hypothesis formulated in [551 HO], despite the fact that they originate 
from the double logarithmic behavior, the singularities at M = j — 2 = — r (r = 2,3, ...) 
may be predicted from the generalized BFKL equation 

(--) + - + \n\) -2*(1) . (6.12) 



-4g 2 V 2/ V 2 

This equation generalizes (16.71) to the case of higher twist operators due to the relation 

\n\ = T — 2. It should be stressed, however, that these higher twist operators are not 
embedded in the sl(2) sector and they still have not been fully identified in the M = 4 
SYM theory. 

Allowing for \n\ to become negative \n\ = — r + l,r = 2,4, ... such that 

uj = M + r-^0 and \n\ + r - 1 = d(r) to 2 + O{to 3 ) , (6.13) 
one can, after replacing the argument of the second psi function by 

1 + l + M ^l + ^ + \n\+to (6.14) 

as explained in [53, [70], derive the following relation from fl 6 . X 2 1) 

7(2^ + 7) = -16g 2 . (6.15) 

Physically, this corresponds to taking into account the double logarithmic contributions to 
scattering amplitudes ~ (aln 2 s) n s~ r+2 . 

Equation (16.151) is a polynomial equation of the second order. One of the solutions is 
unphysical since the anomalous dimension must vanish when g — > and thus 



16g 

7 = -u + u^/i- — 



ur 

(-An 2 ) (-4o 2 ) 2 (-4o 2 ) 3 (-4o 2 ) 4 
= 2 l g 1 -2 { y J +4 l \ ] -10 1 y (6.16) 

io to 6 to b to' 

6.2 Dynamical Test of the Asymptotic Bethe Equa- 
tions 

The predictions (I6.9P and (I6.16P can be compared with the analytic continuations of the 
one-, two- and three-loop corrections (15.221) . (I5.32p . (15.331) and most importantly with the 
analytic continuation the four- loop result. For M = — 1 + to one finds 

y» = 2 (=*£)_ o (=*£)' + (z^V - 2 ti|31±... . (6,7) 
\ to J \ to J \ to J to' 

Around the negative even values M = —2 + to, —4 + to, . . . one derives 

(-V) 2 (~4ff 2 ) 2 | ^-Ag 2 ? 1Q (-4^ 4 
to to 3 to 5 to 



ABA = 2 V ^^ _ 2 V +4 V _ 1Q V 1£_^ (QJg) 

•3 ,,,5 ,,,7 
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One thus infers that the one-, two- and three-loop leading singularities, as derived from 
the asymptotic sl(2) Bethe equations fl5.10p . coincide with the LO BFKL and the double 
logarithmic predictions. On the other hand, the four-loop correction violates strongly the 
pole structure predicted by the BFKL equation. After analytic continuation some of the 
harmonic sums of the seventh order exhibit poles at M — — 1 of the order higher then 
four. According to BFKL equation, however, these poles should cancel each other. This 
contradiction proves unambiguously that the four-loop result is incorrect and consequently 
that the asymptotic Bethe ansatz (15.101) must fail for twist-two operators at four-loop order. 
Moreover, this also confirms that the wrapping interactions must be taken into account at 
this order. 

It turns out that this maximal violation of the BFKL prediction can be easily traced 
back in table 15.11 since after analytic continuation only the first two sums exhibit poles of 
seventh order 



with r = —1, —2, . . .. The expansion (16. 9ft necessitates that coefficients in front of the both 
sums should be equal, while consistency with ( 16. 161) fixes these coefficients to 5. On the 
other hand, the coefficients following from the ABA are 4 and 6 respectively (see table 



6.3 The NLO BFKL Equation and the Generalized 
Double Logarithmic Constraints 



In this subsection we will briefly discuss refinements of the BFKL and the double loga- 
rithmic predictions. These provide additional constraints on the form of the perturbative 
anomalous dimension of twist-two operators and may therefore be used in order to check 
the validity of any yet-to-be-found system of spectral equations of the Af = 4 gauge theory. 

The leading order BFKL equation discussed in subsection 16.11 has recently been gener- 
alized by including sub-leading effects in the Regge kinematics. This was first calculated 
for the QCD in [71] and later the M = 4 part of the QCD answer was extracted in [551170] . 
Using the dimensional reduction scheme the NLO BFKL equation may be written as 




(6.19) 



LO 



(6.20) 



where 




(6.21) 



5(j) 



4x"(7) + 6 C(3) + 2 C(2) X {l) + 4 X {l) x'{l) 




(6.22) 
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The function $(7) is given by the following expression 



* (7)= f JziL 
1 £j(*+7)' 



#(Jfe + 7 + !)-#(!) 



(6.23) 



The perturbative pole structure of the anomalous dimension may be found similarly as 
in the previous subsection. Since the NLO BFKL equation contains information about the 
next-to-leading effects, it predicts the residue of the next-to-leading pole at each order of 
perturbation theory. Thus, expression (16.91) is refined to 

7 = (2 + 0u)f^f\ -(0 + Ow) {^f) +(0+C(3M {^fj (6-24) 
5^ \ 



4C(3) + i C(4).j ^) ±. .. 

This expansion is reproduced through the one-, two- and three-loop results (I5.22p . (I5.32p 
and (I5.33p . At the four-loop order one obtains from (16.241) all together four constraints. 
The first three leading poles must be absent (though, exceptionally, the vanishing of the 
seventh-order pole implies the vanishing of the pole of the sixth-order) and the residues of 
two further poles must coincide with (16.241) . 

In the publication [72] a conjecture has been put forward, according to which the 
inclusion of the NLO and the NNLO corrections to the doube lograithmic scaling amounts 
to modifying equation (I6.15P in the following way 

1 (2uj + 1 ) = -lQ g 2 (l-S 1 uj-(S 2 + ( 2 )uj 2 )-64g\S 2 + ( 2 -S 2 1 ) 

-Ag 2 (S 2 + S^ 2 ) 1 2 . (6.25) 

Harmonic sums in this formula are all functions of (r — 1). It is straightforward to find the 
corresponding perturbative pole structure 
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2 {-4g 2 Y , (6.26) 



with the coefficients q given by 

ci(u) = --S 1 -cu(( 2 + S 2 ) + ..., 

n ^ 1 , 2 Si ( 2 + S 2 
c 2 (uj) = H 5-H h..., 

UJ 6 UJ 1 UJ 

( s 2 6 Si , -4(( 2 + S 2 )+4S 2 + (S 2 + S- 2 ) 
caH = ^5~7JT + ~ 3 + ■••' 

r(A ' 5 , 20 14 (( 2 + S 2 )- 24 Sl-A(S 2 + S- 2 ) 

c 4 {uj) = ^ + — — + + . . . . (6.27) 

UJ 1 UJ° UJ° 
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It should be stressed that the expression (I6.26j) together with (I6.27P determines residues 
of the poles at each negative even value of M = —2, —4, . . ., thus allowing to compare the 
residue functions and not single numbers only. 

Yet another possibility of testing the validity of the spectral equations offer the reci- 
procity relations [73]- [75], which originate from some curious hidden symmetry of the 
anomalous dimension of twist-two operators. Surprisingly, the same relations may be 
found from the string theory side [76] suggesting their non-perturbative validity. 

Very recently, the wrapping correction to the four- loop result, cf. table 15. 1[ has been 
derived in [77] by evaluating the first finite-size corrections to the string sigma model on the 
AdSs x S 5 at weak coupling. This novel procedure was shown to reproduce correctly the 
four-loop anomalous dimension of the Konishi operator calculated in [78] and [79] with the 
usual perturbative methods, and the corrected anomalous dimension of twist-two operators 
has been found to pass positively all tests discussed in this chapter! 
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Chapter 7 



Dynamical Test of the Planar 
AdS/CFT Correspondence 

In this chapter we will discuss interpolating observables, which allow to dynamically verify 
the AdS/CFT duality 

7.1 Interpolating Functions 

The difficulty in proving the AdS/CFT correspondence is mainly caused by the fact that 
the non-perturbative quantization of the M = 4 SYM theory and that of the string theory 
on AdSs x S 5 space is currently not understood. Luckily, in special limits we will discuss, 
one can use the asymptotic Bethe equations to define non-perturbatively gauge theory 
observables, which can be then continued to the results of the string theory obtained from 
the semiclassical quantization. 

One of the simplest dynamical tests of the AdS / CFT correspondence can be performed 
with help of the |-BPS operators, e.g. Tr(^ L ), see section 14. 1.21 These operators are 
primary fields with respect to the superconformal algebra and are annihilated additionally 
by half of the supersymmetry generators Q. It is thus straightforward to derive from (I4.12p 
and (I4.13P the all-loop relation 

A(g)=p = L. (7.1) 

Since the Cartan labels of the su(4) symmetry algebra do not receive quantum corrections, 
it follows from (17. ip that the scaling dimension is protected and non-perturbatively equal 
to its classical counterpart. On the other hand, the string theory on the AdS$ x S 5 space 
has massless excitations which should be identified with the aforementioned BPS operators. 
Such an identification was done in [80J, confirming the validity of the correspondence in 
this simple case. 

The asymptotic integrability discussed in chapterH]is the first step to the non-perturbative 
definition of the N = 4 gauge theory. With help of the asymptotic Bethe equations (I4.60p 
it is possible to calculate perturbatively anomalous dimensions of the operators only up to 
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the order 0(g 2L ). However, for L — > oo one should be able to reproduce the full pertur- 
bative expansion, which could be then continued beyond its convergence radius. Below we 
will discuss this possibility and argue that the results can be compared with the predictions 
of string theory. 

A particularly suitable subsector for this purpose is again the st(2) subsector. As already 
mentioned in ( 5.181) . states of this subsector with the minimal anomalous dimension scale 
for M — ► oo aaj 

A-A = 1 (g) = f^(g) logM + ... . (7.2) 

The universal scaling functional f(g) is defined in the limit L — > oo, L <C logM, and 
thus in the region of validity of the asymptotic Bethe equations. This scaling function 
was defined at the one-loop order in [15J and subsequently investigated at higher-loops 
in [IS] and [H] by means of the asymptotic Bethe equations (I5.10p . It turns out that the 
logarithmic behavior (17.21) is reproduced by the leading magnon fluctuation density and 
that this density can be determined from an integral equation, the BES equation, allowing 
to compute the perturbative expansion of the scaling function to arbitrary order 



fABA{9) = 8 ^-rV + |vrV-16(^vr 6 + 4C(3) 2 ^ 8 

+32 (y^ 8 + lvr 2 C(3) 2 + 40C(3)C(5)) g w ± . . . . (7.4) 

At each order of perturbation theory i only zeta functions or their combinations of the 
order (21 — 1) contribute, moreover all with the overall sign (— 1) £+1 . The convergence 
properties of this oscillating series were studied in [9]. It was shown that the convergence 
radius is equal ■? and that the series admits a natural analytical continuation to the complex 
plane. 

On the string theory side this limit corresponds to the so called "spinning-string" , with 
spin M on the AdS§ space and angular momentum L on S 5 . The classical equations of 
motion in this limit can be exactly solved, as shown in [ST] E2] and [53] El], and the 
corresponding solution can serve as the starting point for the semiclassical quantization. 
The latter was performed to the two-loop order in perturbation theory (in g 1), resulting 
in the following expansion of the scaling function 

3 log 2 K 1 

f{g)=±g 7.5 

7r Att z g 



1 This logarithmic scaling is a special case of the so called Sudakov scaling, see [85] . 

2 It should be stressed that the scaling function of twist-two operators f^ 2 \g) does not need to coincide 
with the scaling function of twist-three operators f^{g), etc. It would be interesting to prove whether 
the universal scaling function f(g) coincides with its finite-length counterparts 

f(g)= ? f i2 \g)= ? f i3 \g)= ? .... (7.3) 

Recently, the above equality was confirmed up to the fourth-loop order for f^ 2 '(g) and f^ 3 \g), see [35], 
[52] and [77]. 
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Here, K = /3(2) denotes the Catalan constant. This expansion may be confronted with 
the strong coupling expansion of the BES equation. We will discuss this further in section 
17^31 

The function f(g) is thus the first interpolating observable of the AdS/CFT correspon- 
dence. It is natural to pose the question whether it is possible to find its generalizations, 
such that asymptotic Bethe equations (I4.60P may once again be used to define it non- 
perturbatively. A detail analysis of the semiclassical quantization of string theory suggests 
that there exists a natural generalization of (17. 2p . see [821 [861 EZ] , when the length is 
assumed to grow logarithmically with M 

M^cx),L-^oo, L = j\ogM. (7.6) 

In the framework of semiclassical quantization g ^> 1 and it is also convenient to introduce 
the following parameter 

4g\ogM Ag 



" L j ■ 

The classical energy of the spinning string in the limit (17. 6p is then given by [82] 



(7.7) 



E = M + L Vl + z 2 + . . . , (7.8) 
and the one-loop correction found in [86] takes the following form 

Ei = ——L= \zVl + z 2 - (1 + 2z 2 ) log \z + Vl + z 2 
vA yl + z 2 L - 1 

-z 2 + 2(1 + z 2 ) log(l + z 2 ) - (1 + 2z 2 ) log Vl + 2z 2 } . (7.9) 

This expression may be also derived from the asymptotic Bethe ansatz at strong coupling, 
see [SB]- Upon identifying the energy with the scaling dimension of twist operators and 
expanding in small j < 1, the following refinement of the scaling behavior (17.21) can be 
found from ([72]) and (173) 



3 1og2 . K 1 j 2 
j 1 

7r " 4n 2 g 8g 



A(g)-M-L = \Ag J*- - j - —- + ±- + . . . ) \ og M +.... (7.10) 



This suggests that the limit (17. 6p could also be defined on the gauge theory side. In sections 
17.31 and 17.41 we will show that this is indeed the case. In the generalized scaling limit (17. 6p 
the anomalous dimension exhibits again logarithmic scaling 

A — M — L = j(g) = f(g,j) logM + . . . , (7.11) 

with the new scaling function f(g,j) depending now on two parameters. The ordinary 
scaling limit (17.21) is recovered for j — » 

f(g,j) =f(g). (7.12) 

In section 17.41 we will derive a closed integral equation for the leading density of the roots 
and argue that the function f(g,j) is analytic in g and j. 
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7.2 Scaling Limit and the BES Equation 



In this section we will discuss extensively the scaling limit ( \7.2\\ . in which the number of 
the excitations diverges M — > oo, while the length L remains finite or grows slower then 
logarithmically with respect to M. 



7.2.1 The Leading One-Loop Density 



It was shown in |49j that the roots in this limit cover densely the interval (-y, y). Since 
the number of the small holes grows slower then logarithmically with M, it follows from 
( 15.60!) ^at they do not form a gap in the magnon root distribution. The integral equation 
for the density (15.691) may thus be written as 

— 

2np (u) + 2n^6(u)-±-^- T -2 [* dv Po{v J =0. (7.13) 

It should be supplemented with the proper normalization condition 

Po(v) = l. (7.14) 



The equation (17.131) was derived in [19] and thoroughly analyzed therein. Contrary to 
naive expectations, even though u(M) ~ 4p diverges for M — > oo, one cannot use Fourier 
transformation for the sake of solving this equation. Rather then that, one defines the 
rescaled density 

Po(u) = M po(u), where u = M u. (7-15) 

This redefinition preserves the measure of the integration dup (u) = dupo(u). Rescaling 
equation (17.131) and using po(u) together with the relations 

' ' n6(v)+o(±=), (7.16) 



' 3 ^(u-^ + ^j^^ + of-^) , (7/17) 



M(u-u') 2 + ^ v ' M{u-u') 2 \M 2 

where V denotes the Cauchy principle value, one derives 

= Att5(u) + 2V [ 2 u' /°( U '\ 9 . (7.18) 
y_i [u — u r 

2 v ' 

It was shown in |1H] that this singular equation is solved by the Korchemsky densit}{§ 

Pk{u) = — log — — — = — arctanh (Vl — 4u 2% ) . (7-19) 



3 The density (|7.19p was first found by G. Korchemsky in [58] while analyzing certain Baxter equation. 
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Since this solution is singular at u = 0, the original density po(u) = -hpo(Mu) should be 
considered as a distribution rather then a function. The one-loop anomalous dimension to 
the leading order in M can be found directly from (15.641) and (17.191) 

7 2 ABA = lf dv^r- = dvp^ + O (M°) = 8 logM + O (M°) . 

(7.20) 

It should be noted that the leading result does not depend on L, which confirms the 
universal scaling behavior at the one-loop order. 

The solution (17.191) can be also partially recovered from (15.521) . Upon rescaling u as in 
(17715]) . and using (l5T68|) together with (l5T3H|) . (15381) and (15321) . one finds 



, s 1 , 1 - 2m 2 . 
p (u) = -log— (7.21) 

Comparing ( 17.21ft with ( 17.191) . one concludes that the result ( 17.211) approximates well the 
exact solution in the interval u £ (— |, |) and only at the boundaries deviates significantly 
from (I7.19p . The discrepancy between (17.211) and (I7.19P is caused by the Imlog term in 
(I5.52p since in the vicinity of u — ^ und u = — \ the approximation (I5.58P ceases to be valid. 
It should be stressed, however, that the difference between (17.191) and (17.211) is negligible 
at the leading order. 

7.2.2 The BES Equation 

To the leading order in M it is possible to recast equation (I5.74p in the form of a linear 
integral equation for the density. For this purpose the large M expansion of the following 
type of integrals needs to be determined 

/■oo 

f(M)= dxh(x) sin (u(M)x) . (7.22) 
Jo 

Here, liniM^oo u(M) = oo and h(x) is a sufficiently smooth function on the interval [0, oo). 
We start by noting that 

lim /(M) = 0, (7.23) 
and therefore the following asymptotic expansion of (I7.22p may be assumed 

oo 

Ca 



f(M) = y . (7.24) 



To determine the first coefficient Co it is enough to integrate (17.221) by parts 

c = \imu(M)f(M)= lim / dxh(x) ( cos (u(M) x) J = h(0) (7.25) 

M—>oc M—>oc Jq y dx J 
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since the boundary terms vanish. Integrating by parts (n + 1) times we find 



n-l 

Co 



c n = Jim u(M) n+1 yf(M) - u (Jy+j j = for eVen n ' ^ ^ 

All odd c n vanish because of the u(M) i— > — u(M) symmetry of (17.221) . 

The derivation of the leading equation is based on the observation that in the scal- 
ing limit the non-linear term in (15.741) may be neglected and equation (15.741) becomes a 
Fredholm integral equation of the second kind. However, not all inhomogeneous terms con- 
tribute to the leading order. Indeed, one can drop all terms except for the one containing 
the fast (universal) holes u^' 2 \ 



7.27) 



it(e*-l) zt(e*-l) 
since the rapidities of the remaining holes vanish in this limit (see ( 15.601) ) 

u p~0 j = 3,...,L. (7.28 
Therefore the equation for the leading scaling function reads 



Z(t) 



47T cos (tu 



z't(e*-l) it(e*-l 

f oo 



-4^-— / dt'e--t'K(2gt,2gt')Z(t'). (7.29) 
e ~ 1 Jo 

Despite the fact that the universal holes receive quantum corrections, these are sub-leading 
at large M and the one-loop asymptotic behavior (15.381) is still valid. One way to motivate 
this is to note that the eigenvalue of the Casimir operator (I5.15P is additively renormalized 
(see [50]) 

1 l(q) 

J ^J = M+-L + ^. (7.30) 

Subtracting the one-loop part from the counting function 

Z(t) = Z (t) + 5Z BES (t) , (7.31) 

and upon identifying, in accordance with ( 15.721) . 5Z B Es{t) with the fluctuation density 
(see pE]) 

5£BE S (t) = 16mg 2 e-> i og(M ) , (7.32) 

one derives 

<W*) = ^TT7 (^(2gt,0)-4g 2 jf dt , K(2gt,2gt')a BES (t , yj , (7.33) 
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where we have used the formula (17.251) to evaluate the integral 

dt'e- t h f K(2gt,2gt')Z (t') = J du{M) J dt' e~? t' K(2gt,2gt') ——Z (t') 

= -iwi K(2gt,0) log M + O(M ). (7.34) 

Equation (17.331) is the celebrated BES equation and was first found in [9]. Let us note 
that the above derivation differs significantly from the original one and confirms that the 
decomposition into the one-loop density and fluctuation density is mathematically well 
defined even non-perturbatively. 

7.2.3 The Strong Coupling Limit 

Directly from the definition (I5.72p and (I7.20p together with (I7.32p it follows that the 
anomalous dimension to the leading order is given by 

7 ABA G?) = (sg 2 -64g 2 jf dt'k{2gt,2gt')a BES {t')j logM + O(M ). (7.35) 

One therefore identifies 

/ ABA G?) = (sg 2 -64g 2 J™ dt'k(2gt,2gt') a BES (t')\ . (7.36) 

Consequently the scaling function f(g) is completely determined by the solution of the 
BES equation. This allows to compare with the string theory result (17. 5p . A systematic 
methoqj of expansion of the BES equation at large values of the coupling constant was 
proposed in [89]. To the first few orders one finds 



r ( 9 ) = 4 9 -^--!Li-... (7.37) 
7T An z g 



in complete agreement with (17.51) . This constitutes one of the first dynamical tests of the 
planar AdS/CFT correspondence. It should be stressed that this would not have been 
possible without the asymptotic integrability. 

7.2.4 Finite Size Corrections 

Beyond the leading order the large M expansion of the anomalous dimension consists of 
the finite part /finite {g, L) and terms that vanish with M — > oo 



7 = /(^)logM + / finite (^,L) + 0^ j — t^tJ . (7.38) 

The somewhat extraordinary order of the first of these terms O ^ log a M j comes from the 
quadratic contribution of the small holes (I5.60|M For the finite part O (M°), however, the 

4 The BES equation at strong coupling was a subject of extensive studies in many different publications, 
see [SD] and references therein. 

5 These corrections are present already at the one- loop level as may be seen by expanding (|5 . 57[) in 
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approximation (17.281) is sufficient and the perturbative expansion of /finite (fi 1 , L) may be 
found directly from (15.741) 



f^(9,L) = (7-(L-2)log2)/0,)-8(7-2L)C(3)^ 
+8 (^A C(3) + (62 - 21 L) C(5)^ £ 6 



. j3-3L)7r 4 C(3) + 5(32-llL)7r 2 C(5) 

15 

+75(127-46L)C(7)J / 
± .... (7.39) 

The function f^Q(g, L) exhibits transcendentality properties similar to the scaling func- 
tion, however in contradistinction to the latter, it explicitly depends on L. Thus, one 
cannot trust the perturbative expansion beyond the wrapping order. Despite this fact 
though, the part strictly proportional to L should not be influenced by the wrapping in- 
teractions suggesting that also f^^(g,L) remains unaffected. 

7.3 The Generalized Scaling Limit 

In this section we will define the generalized scaling limit (17. 6p at the first order of pertur- 
bation theory. 

In this limit the length scales as L = j log M and one concludes from (15.601) that the 
small holes occupy the interval (— c, c). Numerical analysis suggests that the largest of the 
magnon roots is again of the order ±4^ and consequently the Bethe roots condense on the 
interval (— y> ~ a ) U ( a ; if)- ^ e wn ^ P rove later that a = c. 

It is convenient to decompose the corresponding densitj|§ into the singular part Pk(u) 
and the fluctuation density &{u) 

Pm(u) = Pk{u) + d(u) , (7.40) 
where Pk(u) = 1/M px(u/M), see (I7.15p . Upon adding the following term to (I5.69P 

2! dv , PK ^} = 4 logM / ct / + x _ arctan ( M _ q\\ + o(M°) , (7.41) 

J -a {U-Vy + l nM 

one finds that the fluctuation density a{u) must scale as log M/M. This justifies the 
following redefinition 

, , 8 logM , , 
pm{u)=p K {u) — — a(u). (7.42) 



3 We will denote the density in this limit by p m (u). 
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The equation for the density (15.691) may now be transformed into an equation for the 
fluctuation density 

1 j 1 

2 7r cr(tt) (arctan(w + a) — arctan(w — a)) H r- (7.43) 

2ir 8 u + t 

Vi-oo 7a 7 (« - f) 2 + 1 

This should be supplemented with the normalization condition for the density 

du Pm (u) = l, (7.44) 



oo J a 



from which one derives a 

j = — — 8 [ dua{u). (7.45) 
7T J_ a 

Intuitively, the fluctuation density describes perturbations around the Korchemsky density, 
caused by the gap in the root distribution. The size of the gap 2a may be determined in 
the following way. One substitutes (17.451) in (17.431) and subsequently solves the equation 
for a(u, a). Putting the solution into (I7.45P results in the relation j = j(a), which must be 
then inverted. 

The one-loop anomalous dimension may be found immediately using (I7.19P and (I7.42p 

7 2 ABA (j, M) = [ 8 - — arctan (2a(j)) - 16 ( / * ' du + f° du ) I lo g M 

+O(M ) . (7.46) 
This proves that the validity of (17.61) at the one-loop order. 

7.3.1 Fluctuation Density in the Fourier Space 

Equation (17.431) becomes particularly simple in the Fourier space. We rewrite the equation 
in the form 

1 j 1 

a(u) = — - (arctan(w + a) — arctan(w — a)) T (7.47) 

An 2 16 7i u -\- h 



+ 



dv a(v) f a dv a(v) 



7T 1 + (u — v) 2 J_ a it 1 + (it — v) 



Sticking to the conventions in [49], we define the Fourier transformation of the density to 
b<0 

roc 

a(t) = e -3 / due- itu a(u). (7.48) 



7 Every other quantity in the Fourier space is defined through (|5.73[) . 
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It then follows from (17.471) that a(t) satisfies the following integral equation 

a(t) = -J—^ (k h (t, 0; a) - A _ 4 J™ dt' K h (t, t'; a) a(t')) , (7.49) 
with the integral kernel K^(t,t r ] a) defined through 



a 

du cos(tw) cos(t'u) 

-a 



1 _t tcos(at') sin(at) — t'cos(at) sin(at') l_ , . 

e 2 — e 2 . (r.oUj 



2 7T t t 2 -t 
The normalization condition (17.451) is equivalent to 

4a 16 f°° , A . . t sin at .„ 
j = / dta(t)e 2 ^. (7.51) 

71 71 Jq t 

The significant difference between (I7.43P and (17.491) is the fact that the domain of the 
integration of the latter does not depend on the boundary parameter a and the equation 
may be iterated more easily. Once (I7.49P has been solved to the desired order, one can 
find the corresponding anomalous dimension from 



7 2 ABA (J,M) 
logM 



2 

1 arctan 2a (7-52) 

7T 

- 4 dt (&(t) - 4 1 dt' K h (t, t'- a) a(t') J + O (j^j) ■ 



7.3.2 The Density For the Holes 

The above integral equation for the density in the generalized scaling limit was derived 
using magnon roots. However, as pointed out in the section 15.3.21 the system may also be 
equivalently described by holes. Below we will take this point of view and show that the 
resulting equation for the hole density is equivalent to (I7.49p . 

As discussed above, in the limit (17. 6p the small hole roots cover densely the interval 
(— c, c). Similarly to the relation (15.681) . the hole density can be be defined thorough! 

^^Z{u) = 2irp h {u)+o(jj, ue(-c,c). (7.53) 

It is also very natural to normalize the density to one 

dup h (u) = l. (7.54) 



8 In what follows by hole density we understand the density of the small holes u^\j = 3, . . . , L . 



73 



Approximating the non-linear term in (15.521) as in (I5.58j) we find 
+l21og2- — i^{]. + iu)+^{\-iu)) 

+ / ^ (ip(i(u-v))+iP(-i(u-v)))p h (v). (7.55) 



— c 



Moreover, using (I5.38P one gets to the leading order 

2 11 1 f c dv 

PM = (V^+^H^K 2 + J 2n ^ ( u ~ v )))Ph(v) ■ (7.56) 

Upon solving this equation, the one-loop generalized scaling function may be found from 

flESZD 

-ABA/ • M) f c 11 

= jjm, 12 ]og U M ' = 8 + 2 j J d up h (u) ^(-+iu)+^(--iu)-2^(l)) . (7.57) 

Even though (17.561) together with (17.541) constitute the sought-for system of equations, it 
is more convenient for the computational purposes to transform them in a way such that 
the boundary parameter c does not specify the integration domain. Hence, we rescale the 
variable u and the density itself as in (17. 15ft . 

u 

« = Ph{u) = jcp h (u) , (7.58) 

and define the non-singular kernel 

K(u, v) = — ( tp(ic(u — v)) + c {u — v)) —ip( — h icu) — ip( icu)) . (7.59) 

2n \ 2 2/ 

Equation (I7.56P may now be rewritten as 

2 f 1 

Ph(u) = -c+ / dv K(u,v) p h (v) . (7.60) 

K J-i 

One notes that the dependence on the boundary parameter c is hidden in the integral 
kernel. It follows from (I7.54p and (17.581) that the normalization condition for the rescaled 
density is given by the following expression 

j = J dvp h (u). (7.61) 
Finally, the one-loop relation (I7.57P may be rewritten in the following form 

=8 + 2 J dup\{u) + icu) + - icu) -2^(1)). (7.62) 
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The integral equation (17.601) is particularly suitable for the iteration. To the first few orders 
one finds from (17.611) and (17.621) the following expansion of the one-loop generalized scaling 
function 

f(i)U)= = 8-8jlog2 + ^j 3 7r 2 C(3)-^j 4 7r 2 log2C(3) 

+2/ Q7r 2 log 2 2C(3)-^7r 4 C(5)) +0(f). (7.63) 

It turns out that the equation (17.601) can be also very effectively analyzed numerically. In 
figure EH we present the convergence function of the series /(i)(j) = YlkLo fi>k j k > 

r(k) = (/ lfJfc )-i , (7.64) 

for the first six hundred terms. One concludes from this numerical analysis that the 
convergence radius 

r = limsupr(A;) ~ 0.4 (7.65) 

fc^oo 

exists and thus the function (I7.63P is an analytic function in the region \j\ < r. 

It should be noted that due to (I5.68P and (I7.42p the hole density should be considered 
as an analytic continuation of the fluctuation density ct(u) 

2 

j ph(u) = 8a(u) ue(-c,c). (7.66) 

7T 

Using (17.481) one can rewrite the above relation as 

2 g poo ^ 

j ph(u) = / dt a{t) e2 costu . (7.67) 

This relation may in turn be used to prove the equivalence of (I7.49p and (I7.56p . Multiplying 
(I7.49p with ea costu and integrating in t over the positive real axis one finds, using (I7.50P 
together with (17.511) and (17.671) . the equation (17.561) . Thus, we conclude that a = c must 
hold and that there is no gap between the hole and the magnon distributions. 

7.4 The Generalized Scaling Function To All Orders 

In this section we will generalize equation (I7.43P to all orders in perturbation theory, thus 
proving the existence of the scaling function ( IT. 1 If) . 

7.4.1 Derivation 

The asymptotic Bethe equations (I5.10p in the novel limit (17.61) are particularly easy to 
study with the formalism developed in section 15.41 The non-linear term in (15. T4[) may 
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Figure 7.1: The convergence function r(k) for k = 1, . . . , 600 . 



be again neglected. On the other hand, the small holes cannot be disregarded anymore. 
This is due to the fact that their overall contribution is proportional to L = j log M and 
therefore must be taken into account. Thus, equation (I5.74p in the limit (17.61) takes the 
following form 



Z(t) 



2 7T Lea 



it(e* - 1) 



-J (2gt)-Yl 



L 2 7T COS ( t U 



itle* - 1 



L-2 



-4^ 



e.2 



dt'e'^t' dK(2gt,2gt')Z(t') 



e* - 1 jo 

In similarity to (17.3 ip we subtract the one-loop part 

Z(t) = Z (t) + SZ(t) , 
and identify SZ(t) with the fluctuation density a(t) ^ 

5Z(t) = lQTriei ^& logM. 



(7.68) 



(7.69) 



(7.70) 



9 Thc fluctuation density in the M-space cr(u) for u <E (—a, a) must be considered as an analytic contin- 
uation. 
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This allows to rewrite the equation (I7.68P in the form 



e* - 1 



g 2 K(2gt,0)- 



3 J (2gt) 



+ 



1 L 

— Y 



' t l 2 cos(tu^ ] ) 



2 1 L> 

2 logM ^ 

& 3=3 



81ogM 

dt' K(2gt,2gt') e' 1 ' 12 cos(t'<0 



J=3 



(7.71) 



-4<? 2 / dl!k{2gt,2gl!)cr{l!) 



It follows from ( 15. 80ft that the anomalous dimension to the leading order in M can be 
written as 



(1 L f c 



dt 



Ji(2gt) t/2 rti>ii( + ii i.yi. 



2^ 



cos(twj. 



(7.72) 



o 



2#t 



As discussed in the section I5.4[ the equations for the holes to all-loop order (but for 
perturbative values of g) are given by 

Z(u{) = tt(2< + 5-1). (7.73) 

Upon performing Fourier transformation of this equation, one finds the relation 

i 

71 

which after substituting (I7.69P in (I7.73P can be written in the following form 



sin 



(tuQ Z(t)=Tr(2n{ + 5-l) 



(7.74) 



2nn 



4F(uJ 1 fe) ,4 1) ) - 16 lo S M 



dt 



a(t) 



e 1 ' 2 sin(t< J ) , 



t 



(fc)^ 



(7.75) 



with 



costa; 



e 2 — cos t y 
e* - 1 



ip(i(x — y)) + ^ (— « (^ — y))+ip (i (x + y)) + ip (— i (x + ?/)) 



-2^ 



z a; 



(7.76) 



For L = j log M — > oo the hole roots, as shown before, occupy a finite interval around 
zero. Therefore one easily derives from (17.751) the following integral relational for the 



10 The magnon density is related to a(t) through 

2 1 

7T M 



n 2 1 (1) 8 log A/ 



dta{t) e t/2 cos(*u). 
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all- loop hole density ph(u) 

9 q poo 

jp h {u) = ——F\u,u^)-- / dta(t)e^cos(tu). (7.77) 

7T log M 7T Jq 

The first term in this equation, | j^ffa,!!^), is in the limit M — > oo identical with 
(I7.2ip when expressed in the original variable u = M u. Therefore it follows from (I7.19P 
that 

F'(u,v^) = log M + O(M ) ue(-a,a). (7.78) 

The sums in (17.711) to the leading order may be replaced by integrals over the hole density. 
Using (I7.77P and (17.781) one finds the equation 



a(t) 



e* - 1 



-j- J (2gt) + K h (t,0;a)-A [ dt' K h (t,t'; a)a(t') 
8t Jo 

POO 

+g 2 K(2gt,0)-4g 2 / dt' K{2gt,2 g t') a{t') 
Jo 



poo / poo > 

- 4 g 2 / dt' t' K(2gt, 2gt') K h (t', 0; a) - 4 / dt" K h (t', t") a(t") 
Jo V Jo / 

(7.79) 

with Kh(t,t'; a) defined in (17.501) . From the normalization condition for the hole density 

dup h (u) = l (7.80) 



one derives the same relation as in the one-loop case (I7.5ip 

4a 16 f°° , „ , s t sin at ,„„„\ 

j = / dta(t)e*——. (7.81) 

n n J t 

Finally, putting (17.811) in ( 17.791) one finds the desired integral equation for the leading 
density in the generalized scaling limit (17.61) 

a(t) = (fc(t, 0) - 4 jf dt'JC(t, t') a(t')j . (7.82) 

The new kernel JC(t,t') is an intricate combination of the hole an magnon kernels 

* - . „ * . , . Jn(2qt) sin at' *' 

JC(t,t') = g 2 K{2gt,2gt')+K h (t,t'-a)--^^——e2 

t Zttv 

poo 

-Ag 2 I dt"t" K(2gt,2gt") K h (t",t';a) . (7.83) 
Jo 

It is interesting to note that equation (I7.82p takes structurally similar form to the BES 
equation (17.331) . The leading anomalous dimension (I7.72p may be expressed through a(t) 
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with help of (I7.77p . which in conjunction with (17. lip leads to the following expression 



f AaA (g,j) = 8 g 

- 8 



oo 

dt 



o 2gt 

J x {2gt) 



a(t)-4t I dt' K h (t,t';a(j))a(t') 
o 



16 <r(0) + 



J_ 
16 



(7.84) 



Since the wrapping interactions may be neglected in the limit in question, the above formula 
proves the existence of the generalized scaling function (17.111) to all orders in perturbation 
theory. This novel function depends on two variables and thus it is possible to define two 
infinite families of scaling functions by expanding in g or j respectively 



f(g,j) = f(g) + J2f (n) ^)j n = /d)(i) + Y,hn)U)g 



2n-2 



(7.85) 



n=l 



n=2 



In particular, the function f^(g) coincides with (17. 39ft after setting L = j logM and 
keeping the leading terms in the expansion only 



f (1 \g) 



-f(g) log 2 + 16s 4 C(3) -^-vr 2 C(3) + 168 C(5) 
+S 8 (| vr 4 C(3) + f vr 2 C(5) + 1840 C(7)) + . . . . 



(7.86) 



This suggests that the L dependent part of (I7.39P is not influenced by the wrapping ef- 
fects. Another curious observation is the vanishing of the f^ 2 \g) function to all orders in 
perturbation theory. An indication of such behavior on the string theory side was observed 
in [91]. 

Interestingly, for j g it was conjectured in [ST] that the generalized scaling function 
(17. lip can be computed using the results of the 0(6) sigma model. Recently, it was 
confirmed [92] that the equation (17.821) together with (17.811) are equivalent to the TBA 
equations of the non-linear 0(6) sigma model. Please refer to |93j-[95] for the further 
development on this subject. 
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Chapter 8 

The Hubbard Model 



In this chapter we will explain an interesting relationship between the su(2) asymptotic 
Bethe equations without the dressing factor and the spectral equations of a well-defined 
short-range integrable model, the Hubbard model. 

8.1 The 5u(2) Sector and the Hubbard Model 

The su(2) sector is one of the simplest dynamical sectors in the M = 4 gauge theory. The 
operators of this sector are composed of two scalar fields only, which we will denote by Z 
and X. Conventionally, we choose the reference vacuum to be 

\Z L ) = TiZ L . (8.1) 

Excited states are obtained by putting X fields into the trace 

Tr (X M Z L - M ) + . . . , (8.2) 

followed by the diagonalization of the su(2) dilatation operator, which is a consistent 
truncation of the complete dilatation operator to this sector. The Cartan weights of the 
primary state corresponding to the diagram 14.41 are given by 

{A , s u s 2 , q uP , q 2 , B, L} = {L, 0, 0,M,L- 2M, M, 0, L} , (8.3) 

from which, according to (14.671) . the following excitation pattern emerges 

{K u K 2 , K 3 , K A , K 5 , K 6 , K 7 } = {0, 0, 0, M, 0, 0, 0} . (8.4) 

The corresponding all-loop asymptotic Bethe equations^ take particularly simple form 

^ = II k Uj + l ew(2ie(u k , Uj )). (8.5) 

1 In the literature these equations are known as the BDS equations 40J. 
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At the one-loop order the above equations reduce to the spectral equations of the 
Heisenberg spin chain (I3.52p and consequently the dilatation operator of su(2) may be 
identified with the Hamiltonian of the XXXi spin chain (13.421) . The higher corrections 
to the dilatation operator were studied in [39], [40] and [IT], and under some assumptions 
were derived up to the five-loop order. It turns out that the £-loop dilatation operator acts 
simultaneously on the (£+ 1) neighboring lattice sites. Moreover, these corrections become 
very complicated beyond the first few orders. Therefore it seems hopeless to guess its all- 
loop form solely from the perturbative expansion. In order to reveal the hidden connection 
to a different integrable model, we will investigate in detail the antiferromagnetic state of 
the su(2) sector. 

We start by noting that the S-matrix on the right hand side of (|8.5|) . upon neglecting 
the dressing factor, depends only on the difference of the rapidities. Since the dressing 
factor starts to contribute from the fourth-loop order, we will neglect it in what follows 
and discuss this approximation at the end of this section. The antiferromagnetic state is 
a maximally filled state. In the case of compact spin chains like XXXi this implies that 

M mSLX = § @- The energy of this state in the case when the length of the chain becomes 
thermodynamically large was computed for the XXXi spin chain as early as in 1938 by 
Hulthen [96J. He observed that the corresponding Bethe equations in the limit L — > oo 
can be written in the form of an integral equation, which in turn may be solved by Fourier 
transformation. Below we will apply this method to the deformed equations (18.51) . 

Taking the logarithm of the equations (18.51) and introducing the root denisty, one finds 

f°° , , p(u') d , x+(u) , , 

2 ^W + 2 //" (M -V + i = ^ log FR (8 - 6) 

since in the Arctan branch of the logarithm the mode numbers are uniformly distributed. 
It follows directly from (I4.59P that 

i— log — ^4 = —!= ■ 8.7 

du x-(w) + -Ag 2 y/(u - i/2f - Ag 2 

The integral kernel in (18.61) depends on the difference of the variables only, and the equation 
may thus be solved by Fourier transformation 

[°° dt cos (to) M2gt) 

" {U) = L aT cosh (I) ■ (8 ' 8) 

The corresponding anomalous dimension may be found immediately from (14.681) 

i{g) 



. dupiu) — — — + 0(L 

2g*L ^ ; \x+{u) x~{v> 1 



oo 



4/ ±MW)M2 g t) 



2 Equations (|8.5p . under the assumption that j« 1, may be considered as an all- loop deformation of 
the Heisenberg spin chain. 
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It is now very interesting to note that equations (18.81) and (18.91) are well-known results 
of the solid state physics. They describe respectively the density and the energy of the 
antiferromagnetic state of the Hubbard model. The integrability of this model was proven 
by Lieb and Wu in [97], where also the corresponding Bethe equations were derived. 

The Hubbard model is a dynamical, short-range model of N electrons on L lattice sites. 
Due to the Pauli's exclusion principle, there are four possible states on each lattice site 

1. no particles, 

2. spin-up electron f, 

3. spin down electron j. and 

4. double occupied state with spin-up and spin-down electrons J:=TI- 

In what follows, we will consider the half-filled case N = L. The Hamiltonian of the 
Hubbard model consists of the kinetic part that forces the electrons to jump between 
different sites and the potential part, which according to the value of U corresponds to 
repulsive or attractive force 

L L 

^Hubbard = -t^2^2 ( C L C m> CT + 4+l,* C i,<r) + t U Yl 4,^4,1^,1 ■ ( 8 ' 10 ) 
t=l a=1,l i=l 

The operators c\ a and Cj CT are canonical Fermi operators obeying 

{ci,„,c jiT } = {4 CT , c t T } = 0, (8.11) 
{c£,<7,ct T } = S i:j S aT . (8.12) 

We assume the system to be closed and thus we identify 

c L+l,a = Cl,a, C L+1,(T = C l,a (8.13) 

for a =|, \. The Hamiltonian is invariant with respect to the su(2) transformations 

[^Hubbard, S a ] =0 a = +,-«, (8.14) 

with S a = J2i=i Si 1 - This allows to classify the spectrum according to the eigenvalues of 
the total spin an its z component. 



The precise correspondence between (18. 8p . (18.91) and the corresponding results of the 
Hubbard model is established under following identification of the parameters 

t = ~, U = -. (8.15) 

2^ g 

This, however, suggests that perhaps all states of the BDS spin chain are up to the wrapping 
order also eigenstates of the Hubbard model Hamiltonian. It turns out that in the case of 
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odd values of L the BDS equations indeed diagonalize the Hamiltonian (18. 10p . while for the 
even values the kinetic terms of the Hamiltonian need to be multiplied by additional phase 
factors. This is equivalent to coupling the Hubbard model to a homogeneous magnetic 
field and the aforementioned phase factors may be then considered as the Aharonov-Bohm 
phases 

1 L 1 L 

^ = ^EE c L c i+v + <W**) - E c WsW > ( 8 - 16 ) 

- J j=l a=tl y 3=1 



with 



T = i = 0, 0=^ mod2|. (8.17) 

The perturbative expansion at small values of g corresponds accordingly to the strongly 
coupled (U 3> 1) Hubbard model. It was shown in [98J that the perturbative expansion of 
the Hamiltonian (I8.16P at large values of U = 1/g coincides up to three orders in g 2 with 
the perturbative expansion of the dilatation operator found in [40J . We will not repeat this 
calculation here, but instead show in the next section that the BDS equations (18.51) can be 
derived in the asymptotic region from the Lieb-Wu equations. 

The vacuum state of the Hubbard model is annihilated by the c^, c^j operators and is 
a tensor product of L empty lattice sites 

\0) L = |0)® |0)®...|0> . (8.18) 

" . ' 

L 

The elementary excitations on this vacuum are the spin-up and spin-down electrons. On the 
other hand, the double occupation must be considered as a composition of the elementary 
excitations. The interaction between the constituents of such a composed state is according 
to (18.161) repulsive. Since we confine ourselves to the half-filled case, it is more convenient 
to choose the BPS vacuum as the reference vacuum 

\z L ) = | TT • • • TT) = AA ■ ■ ■ ■ (8-19) 

It is easy to show that this state is also a zero energy state of (I8.16p . The disadvantageous 
feature of the vacuum (I8.19p . that is the fact that it is not annihilated by the annihilation 
operators, may be easily overcome by performing a particle-hole transformation 

^=> t, (8.20) 

1 I ■ (8.21) 

After this transformation the spin-up electrons are considered to be empty lattice sites, 
while the empty lattice sites become double occupied states. Explicitly, we identify 

c j,o — A,t ' c l° = c i'T ' (8.22) 
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The algebraic relations between Cj ,Cj j0 ,cl + and c^j follow directly from the anitcommu- 
tation relations ( 18.1 ip and ( 18.121) . Accordingly, the dual Hamiltonian can be written as 

L L 

^ = ^ £ £ ( e * <Wv + <W*>) - 0^2 £a - cWjc],^,! , (8.24) 

9 j=l £r=o,I 9 i=l 

with the phases 0j = 0j_ and O = 7r — 0^ being different from the one for the origi- 
nal exciations. Comparing (I8.16P and (18.241) one finds that performing the particle-hole 
transformation amounts to the substitution 

M 

H(g; 0, 0) - -H(-g; vr - 0, 0) - — , (8.25) 

where M denotes the number of the double occupancies. The sign in front of the potential 
part of the Hamiltonian (I8.24p is opposite to the one in (I8.16p . which means that the 
holes and the double occupancies attract each other and form bound states, the spin down 
electrons. 



8.2 The Lieb-Wu Equations 

The Hamiltonian (I8.10p . as shown in [97], is integrable and can be diagonalized with the 
help of the Bethe ansatz. The same holds for Hamiltonians with arbitrary Aharonov- 
Bohm phases, which were studied in [HS]. In the case of half-filling the Bethe equations 
corresponding to (I8.16P are given by (see |99j) 

e ., J = n %-2 g smfe^)-*/2 
J-J^ Uj - 2g sm(q n + 0) + i/2 

*u k -2gsm(q n + <f>)+i/2 ^ u k -u s + i ^ 
M u k -2gsm{q n + (p) - % 2 . | A u k -Uj—i 

Here, the parameter M denotes the number of the spin-down electrons. The energy of a 
state is a function of all q n 

1 L 

E= - Vcos(g n + 0) . (8.28) 

y n=l 

In the limit g — > equations ( 18.261) and ( 18. 27ft become the momentum constraint (after 
setting e^ nL = 1) and the Bethe equations of the su(2) spin chain respectively. The 1/g 
correction to the energy ( 18. 28ft is zero due to 

e nnL = 1 ^ q n = ?l( n -i) + o(g), n = l,...,L. (8.29) 
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To obtain the O(g ) term one must determine the 0(g) corrections to the roots q n . Since 
the one-loop roots are known, this is merely a linear problem. Upon solving (18.261) to the 
order 0(g), one finds 

- 1 

E = J2 T7T + °(9) , (8-30) 

k=l U k + 4 

which is the expression for the energy of the XXX spin chain, cf. ( 14.41 [) . Knowing the 
one-loop solutions q n and u k , one can solve equations (I8.26P and (18.271) perturbatively. In 
general, for 1 <C L < oo, it is very complicated to solve the Bethe equations exactly. On 
the other hand, the BDS equations are strictly speaking only valid for L — > oo. Therefore 
it is convenient to dualize equation (I8.26p . Introducing x = e^™ one concludes that the 
equation (I8.26P in its polynomial form has in total L + 2M solutions. Using the remaining 
2M roots, one can rewrite (18.261) and (I8.27P as 



eiqnL ^ u j -2gsHc ln -4>)-i/2 



. ;l Uj - 2gsin(q n - 0) + i/2 



' u k - 2gsm(q n -<!>) + i/2_ rr u k - Uj + i fc = 
n=[ u k ~ 2 9 sm(g„ -(f))- i/2 a 4, u k -Uj-i 



The energy in terms of the dual roots can be found to be 

2M 



E 



M 



V 9 „ , 



- J^cos(g n - 0) . (8.33) 



Comparing (I8.33j) to (I8.28P together with (I8.25P one infers that the dual solutions q n ,n = 
1, . . . 2M diagonalize the dual Hamiltonian. 

The advantage of the equations (18.31 j) and (18.321) . as opposed to (I8.26P and (I8.27p . is the 
independence of the former on the length L. Since the elementary dual excitations attract 
each other, it is natural to assume that bound states (|= o + J) will be energetically 
favored. Usually, the bound states manifest theirselves in the exponential localization 
of the wave function. Therefore, for the lowest states in the energy band we make the 
following ansatz for the 2M dual roots^l 



q n -4> = -sgn(p n ) + y + i(3 n , (8.34) 
Qn+M -<P = -sgn(p n ) + — - - i/3 n , j3 n > , n=l,...,M, 

3 The separation of the -|sgii(p„) factor is made in order to simplify the calculations. The variable p n 
will correspond to the momentum of a magnon p n £ (— tt, it) so that (q n — <f>) S ^). Interestingly, 
there exist no solutions of (|8.3ip and (|8.32[) with (q n — (f>) S (— tt/2, 7r/2). Such solutions would correspond 
to magnons with negative energies. 
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with p n e M. denoting the bound state momentum. Substituting (18.341) and (I8.35P into 
(18.3ip one finds that when L> 1 the left hand side of this equation vanishes like e~@ nL for 
n — 1, . . . , M and diverges like e^ nL for n — M + 1, . . . , 2M. This is true independently of 
the value of M, and therefore the following relations must hold up to exponential corrections 

u n - i/2 = 2g sm(q n - 0) + 0(e~^ L ) , u n + i/2 = 2g sm(q n+M - <j>) + 0(e~^ L ) . (8.35) 

Putting ( 18.351) into the equation (18.321) one finds that the latter is trivially satisfied. The 
both conditions (18.341) and (18.351) may be compactly written as 

u n ± i/2 = 2g sgn(p„) cos(^ =F ifi n ) + 0(e~^ L ) . (8.36) 

Solving ( 18.361) for p n and (3 n results in 

sinh/3 n = — =- = - — , 1 P . , (8.37) 

Ag sgn(p n ) sin ^ 4#|sm^| 



u n = 2g sgn(p n ) cos y cosh (3 n = - cot y W 1 + 16^ 2 sin 2 y . (8.38) 



The second relation is equivalent to the all-loop dispersion relation ( 14.601) 



ip 



x + {u) 
x~{u) 



(8.39) 



In order to eliminate q n from ( 18.311) one substitutes the relation ( 18.351) into ( 18.3 If) and 
subsequently multiplies the equation for u n with the equation for u n+ M- The divergencies 
appearing for L — > oo cancel each other and the finite part can be written as 

M _ 

e i{ Pn +2<t>+s n ^)L = _ J-J U n Uj + I ^ (8.40) 
■ -. -i U n Uj % 

Using ( 18.171) and (18.391) one confirms that (18.401) is equivalent to the BDS equations (18.51) . 
The energy (I8.33|) reads 



M 

I ■ 

E 



~ Yl ( C0S (^ -<!>) + C0S (<?n+A/ - #)) - „-2 ( 8 - 41 ) 



M M 



-y^lsin— cosh/5 n - = I — —, — r ; — rl 

g^\ 2 ^ n 2g* ^\x+(u n ) X -(u n )J 



n=l n=l 



M ^( l l A j(g) 

2g 2 



where we have made use of (I8.37j) and (I8.39[) . The above construction can be generalized 
to the case of complex momenta p n . The argument of the sigma function in (18.341) and 
must then be replaced by the real part of p n . 
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For the purpose of derivation of the BDS equation we have assumed that the holes o and 
the double occupancies | form bound states. Though this is true for the low-energy states, 
there exist states for which one or more bound states are split up and the corresponding 
momenta are real. Their presence may be also noted by a simple counting. In the Hubbard 
model there are roughly 4 L /L cyclic states, while in the Heisenberg or BDS spin chain the 
number of cyclic states amounts to 2 L /L. This discrepancy is explained by the existence of 
the excited states. One should note, however, that the energies of these states are according 
to (I8.33P non-perturbative 

E ~ \ , (8.42) 

g 2 

and therefore the corresponding operator cannot be defined in the perturbative gauge 
theory. The accuracy of the equation (I8.40p is related to the accuracy of the solution 
(I8.35p . It follows immediately from (I8.37P that for perturbative values of g 

(3 n ~ - log g + O(g ), (8.43) 

and thus the solutions (I8.35P receive corrections starting from the order e~ l3nL ~ g L . The 
energy (18.411) . on the other hand, due to the expansion 

cos(x + e) = cos(x) — sin(x)e — - cos(x)e 2 + 0(e 3 ) , (8.44) 

and the fact that q n and q n+ M are complex conjugated to each other, leeds to a correct 
result up to the order (D(g 2L ~ 2 ). According to the relation 7(g) = 2g 2 E(g) this corre- 
sponds precisely to the wrapping order. This is a strong indication that there exist an 
integrable short range model, which correctly captures the wrapping interactions and in 
the asymptotic region leeds to the asymptotic Bethe equations (18.51) . It should be noted, 
however, that the usual Hubbard model may not be a candidate for such model since it 
leads to the trivial dressing factor <r(u,v) = 1, cf. (18.401) . This contradicts both the pro- 
posed crossing symmetry [8] (cr(u, v ) = 1 does not satisfy the crossing equation) and the 
explicit perturbative calculation [47j . It is, however, likely that a suitable deformation of 
the Hubbard model would lead to the "dressed" BDS equations. 

The hypothesis that the dilatation operator of the Af = 4 SYM theory is equivalent to a 
Hamiltonian of a short-range integrable model and that the observed long-rangeness is only 
an artifact of the perturbative expansion is very appealing, although several issues need still 
to be understood. Most importantly, it is not known whether such a short-range system 
with both bosonic and fermionic elementary excitations can be found since, in general, 
unlimited number of bosons can occupy each lattice site preventing the factorization of the 
scattering into a sequence of two-body processes. Another complication is the presence of 
the dressing phase, which has a very complicated transcendental structure, see [9]. 

It is also interesting to note, as it was found in |100j . that the S su ( 2 \2) S-matrix may be 
identified with the R-matrix of the Hubbard model. However, this observation does not 
seem to be related with the context in which the Hubbard model was introduced in this 
section. 
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